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INTRODUCTION 

The  mechanism  by  which  heat  or  mass  is  dispersed  in  tubular 
flow  reactors,  and  other  chemical  process  equipment  through  which 
a  fluid  is  flowing,  has  recently  gained  a  great  deal  of  attention,. 
In  order  to  solve  the  problem  of  determining  the  conversion  and 
mixing  characteristics  of  a  tubular  flow  reactor  it  is  necessary 
to  consider  a  rather  complicated  system  of  partial  differential 
equations,,   The  system  is  complicated  because  the  velocity  profile 
is  not  flat  and  because  there  are  not  only  gradients  of  concentra- 
tion but  also  of  temperature  and  pressure 0   To  obtain  the  general 
solution  of  this  problem  analytically  is,  of  course,  out  of  the 
question  and  the  numerical  solution  is  tedious  and  requires  a  very 
large  computer „   At  the  same  time,  however,  there  is  the  larger 
question  of  whether  from  a  practical  point  of  view  the  effects  men- 
tioned above  are  important  enough  to  justify  the  effort  needed  to 
solve  the  problem,,   Therefore,  it  would  be  interesting  to  obtain 
some  estimate  of  the  magnitude  of  some  of  these  effects 0 

In  an  attempt  to  exploit  a  technique  which  might  give  a  partial 
insight  into  the  behavior  of  tubular  flow  reactors,  a  flat  velocity 
profile  and  a  concentration  gradient  due  to  dispersion  only  is  as- 
sumedo  Earlier  investigators  introduced  the  concept  of  longitudi- 
nal diffusion  to  describe  the  phenomenon  of  dispersion  in  the  di- 
rection of  flow  of  the  mean  fluid  stream,  and  the  concept  of  rad- 
ial diffusion  to  describe  dispersion  in  the  direction  perpendicular 
to  the  mean  fluid  flow0   This  work  will  be  limited  to  the  study  of 
axial  dispersion,, 


A  number  of  models  have  been  proposed  for  the  dispersion  phe- 
nomenon associated  with  flow  through  tubular  flow  reactors,  fluid- 
ized  beds  and  packed  beds,  for  example,  the  researches  of  Danckwerts 
(2)  and  Levenspiel  and  Smith  (12)  for  flow  through  tubular  flow  re- 
actors by  the  diffusion- type  model,  the  work  of  Carberry  and  Bretton 
(1)  for  flow  through  fixed  beds  by  the  cell-model,  the  work  of 
Miyauchi  (111)  by  a  cell  with  back-mixing,  and  the  diffusion-dead 
space  model  by  Koump  (9)  and  Turner  (20),  All  of  these  models  ideal- 
ize the  three-dimensional  problem  into  a  one-dimensional  transport 
equation,. 

The  most  commonly  used  model  is  the  diffusion-type  model0  This 
model  assumes  that  the  overall  rate  of  axial  dispersion  is  propor- 
tional to  the  bulk  concentration  gradient  of  the  diffusing  component 
in  the  axial  direction  and  is  measured  by  the  spreading  of  a  concen- 
tration wave  in  an  axial  direction,,  This  model  is  also  used  in  this 
research  to  investigate  the  effect  of  axial  mixing  on  the  perform- 
ance of  isothermal  tubular  flow  reactors.  With  this  model,  a  bound- 
ary value  problem  of  a  second-order  partial  differential  equation  is 
obtained  when  a  material  balance  is  made  with  respect  to  a  traceable 
component  in  the  fluid. 

In  general,  a  change  in  boundary  conditions  may  have  a  signifi- 
cant effect  on  the  behavior  of  the  solution  of  a  partial  differential 
equation,,   Hence,  in  order  to  avoid  the  possibility  of  error  due  to 
the  improper  choice  of  boundary  conditions,  one  must  prove  that  the 
boundary  conditions  used  in  solving  a  particular  initial  or  bound- 
ary value  problem  are  realistic,  or  that  the  solution  of  the  partic- 
ular differential  equation  is  not  sensitive  to  the  choice  of  bound- 
ary conditions. 


Various  sets  of  boundary  conditions  have  been  suggested  for 
use  with  the  d If fusion-type  model  to  predict  dispersion  phenomenon 
associated  with  flow  through  tubular  flow  reactor.   Danekwerts  (2), 
for  example,  introduced  a  set  of  boundary  conditions  which  satisfy 
the  continuity  requirement  at  both  ends  of  a  reactor.   Many  other 
investigators  introduced  sets  of  boundary  conditions  which  do  not 
conserve  the  continuity  requirement  at  the  inlet  of  reactor0   Hul- 
burt  (8),  for  example,  introduced  a  set  of  boundary  conditions  which 
neglect  the  diffusion  at  the  inlet  of  the  reactor,   Levenspiel  and 
Smith  (12)  further  approximated  a  reactor  of  finite  length  by  a  re- 
actor of  infinite  length.   Although  both  of  these  two  sets  of  bound- 
ary conditions  are  convenient  for  estimating  the  relative  importance 
of  reaction  and  diffusion  rates,  there  has  been  no  attempt  made  to 
justify  numerically  the  quantitative  relationships  among  the  solu- 
tlons  obtained  by  use  of  these  sets  of  boundary  conditions .  Another 
possible  set  of  boundary  conditions  which  meet  the  continuity  re- 
quirement at  both  ends  of  the  reactor,  but  hitherto  unknown,  is  in- 
vestigated in  the  course  of  this  research. 

The  present  study  represents  an  effort  to  increase  the  theo- 
retical knowledge  of  isothermal  tubular  flow  reactors  by  investiga- 
ting the  effects  of  axial  dispersion  on  the  performance  of  a  reactor 
based  on  a  diffusion-type  model.   The  differences  and  the  relation- 
ships among  the  steady-state  response  curves  against  the  various 
types  of  concentration  pertubation  obtained  by  use  of  the  different 
sets  of  boundary  conditions  are  numerically  evaluated.   Theoretical 
response  curves  of  exit  concentrations  against  impulse,  step  input 
and  frequency  response  functions  are  presented  in  order  to  illus- 


trate  the  method  of  determining  mean  dispersion  coefficients  from 
available  dynamic  experimental  data0 


LITERATURE  SURVEY 

As  indicated  in  the  previous  section,  different  models  of  fluid 
flow  and  various  sets  of  boundary  conditions  have  in  recent  years 
been  proposed  in  order  to  predict  the  effects  of  axial  mixing  on  the 
performance  of  such  process  equipment  as  packed  beds,  fluidized  beds 
and  tubular  flow  reactors.   Table  1  summarizes  the  work  of  previous 
investigators o   Included  in  the  table  are  the  mathematical  models 
and  boundary  conditions  employed,  the  process  equipment  and  experi- 
mental technique  used,  and  the  method  of  analysis  of  the  data. 

Danckwerts  (2)  presented  a  general  discussion  of  how  distribu- 
tion functions  for  residence  time  can  be  defined  and  measured  for 
pipes,  packed  beds  and  reactors  from  dynamic  step-input  experiments 0 
He  also  reported  the  performance  of  a  steady-state  flow  reactor  based 
on  the  diffusion-type  model  with  a  first-order  chemical  reaction  us- 
ing a  set  of  continuity  preserving  boundary  conditions  at  both  ends 
of  the  reactor0   The  effect  of  axial  mixing  on  the  first-order  con- 
version is  also  compared  with  that  of  a  plug  flow  reactor,.  The 
boundary  conditions  Danckwerts  introduced  will  be  denoted  as  B.C0(I)o 

Hulburt  (8)  investigated  steady-state  first-order  kinetics  us- 
ing a  diffusion  model  to  study  the  effect  of  axial  dispersion  on  the 
conversion,,   He  presented  a  set  of  boundary  conditions  which  neglect 
the  effect  of  diffusion  at  the  reactor  entrance e   This  set  of  bound- 
ary conditions  will  be  denoted  as  B.C. (II). 

An  analysis  of  these  two  boundary  conditions  was  presented  by 
Wehner  and  Wilhelm  (21;)  for  a  steady-state  flow  reactor  with  axial 
diffusion  and  first-order  reaction,,   They  have  indicated  that  the 
Hulburt  solution  (BoCoII)  is  convenient  for  estimating  the  relative 
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Importance  of  reaction  and  dispersion  rates,  but  is  not  consistent 
with  the  limit  of  a  stirred  tank  reactor  with  infinite  diffuslvity 
in  the  reaction  section,,   They  have  investigated  the  validity  of 
these  boundary  conditions  by  simultaneous  solution  of  three  differ- 
ential equations  for  the  reaction  section,  fore  and  after  sections. 
This  investigation  was  based  on  the  fact  that  axial  dispersion  may 
occur  in  the  latter  two  sections 0 

Levenspiel  and  Smith  (12)  provided  an  anlytical  soltuion  to 
the  problem  of  fluid  mixing  In  a  continuous  flow  system  without 
chemical  reaction  using  the  diffusion-type  model  but  with  a  differ- 
ent set  of  boundary  conditions.  They  presented  a  set  of  approximate 
boundary  conditions  by  assuming  that  the  effect  of  diffusion  at  the 
reactor  entrance  is  negligible  so  that  a  discontinuity  exists  at 
this  boundary o   They  further  approximated  a  reactor  with  finite 
length  by  a  reactor  with  Infinite  length.  Based  on  this  set  of 
boundary  conditions  they  showed  two  methods,  a  single  point  method 
(maximum  concentration)  and  a  method  based  on  the  variance  of  the 
concentration  distribution  {£    ),  for  evaluating  the  axial  disper- 
sion coefficient  from  experimental  measurements.   They  also  worked 
out  some  examples  to  determine  the  condltons  of  applicability  of 
the  model.  This  set  of  boundary  conditions  will  be  denoted  as  B,C0 
(III), 

In  a  short  note,  van  der  Laan  (21)  discussed  an  aspect  of  the 
work  reported  by  Levenspiel  and  Smith  (12),   While  Levenspiel  and 
Smith  (12)  treated  the  diffusion-type  approach  by  restricting  the 
solutions  to  infinite  lengths,  van  der  Laan  (21)  treated  the  dif- 
fusional  type  flow  in  a  more  general  manner  so  as  tc  include  the 


case  of  finite  length  and  of  varying  dispersion  by  applying  suit- 
able boundary  conditions 0   He  treated  the  general  solution  by  means 
of  a  Laplace  transformation  and  avoided  the  back  transformation  by 
applying  momentum  analysis,  that  is,  the  method  of  variance  of  the 
concentration  distribution,,   An  illustration  was  given  for  the  eval- 
uation of  the  dispersion  coefficient  using  the  results  of  calcula- 
tions of  mean  residence  time  and  variance,, 

One  of  the  pioneer  investigations  of  axial  dispersion  in  pro- 
cess equipment  was  that  of  Kramers  and  Alberda  (10).   They  investi- 
gated axial  dispersion  for  the  flow  of  water  through  a  bed  of  Raschig 
rings o  A  frequency  response  technique  was  used,  and  the  results 
were  analyzed  on  the  basis  of  the  assumption  that  the  rate  of  axial 
dispersion  can  be  described  by  a  diffusion -type  model  with  B.C. (I), 
In  order  to  avoid  the  complexities  involved  in  obtaining  an  exact 
analytical  solution,  they  have  suggested  that  the  performance  func- 
tion of  the  system  be  approximated  by  n  completely  mixed  vessels 
connected  in  series 0   The  feasibility  of  this  approximation  is  open 
to  question,  except  for  the  system  involving  a  small  degree  of  dis- 
persion,  Kramers,  Honig,  and  Croockweit  (11)  extended  their  work 
to  the  study  of  longitudinal  dispersion  in  liquid  flow  by  using 
B«,Co(II)o   Actual  derivation  of  the  analytical  solution,  however, 
was  carried  out  by  means  of  B.C. (III). 

Yagi  and  Miyauchi  (25)  investigated  the  operational  character- 
istics of  continuous  flow  reactors  in  which  the  reactants  are  mix- 
ing and  described  the  dispersion  phenomenon  in  terms  of  mean  longi- 
tudinal dispersion  applying  the  diffusion- type  model  and  B.C. (I). 
An  impulse  function  was  used,  and  the  results  were  analyzed  by  means 
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of  a  single  point  method  (l0e.s  maximum  points).   The  data  of  Ebach 
(l\.)    on  impulse  experiments  were  analyzed  on  the  same  basis  as  those 
of  Yagi  and  Mlyauchi  (2lj.).   It  is  known  that  the  use  of  a  single  da- 
bum  in  drawing  conclusions  for  an  entire  experiment  may  lead  to  seri- 
ous error  because  of  the  probable  experimental  error,  Although  the 
data  shown  in  the  work  of  Ebach  (if)  gave  close  agreement  between  im- 
pulse experimental  values  and  calculated  values  at  low  degree  of 
dispersion,  the  feasibility  of  using  peak  values  of  their  impulse 
data  at  high  degree  of  dispersion  leaves  much  to  be  desired0 

In  experiments  with  the  frequency  response  technique ,  Ebach 
(k)    used  amplitude  data  alone  to  analyze  the  mean  dispersion  coef- 
flciento   He  stated  that  the  axial  dispersion  coefficient  calcula- 
ted from  phase  shift  data  gave  poor  accuracy  and  reproducibility,, 
The  reasons  for  this  fundamental  fault  were  clearly  illustrated  in 
his  paper*   Strang  and  Geankoplis  (18)  studied  the  axial  dispersion 
of  liquids  in  packed  beds  by  the  same  frequency  response  technique „ 
The  diffusion  model  was  used,  and  the  data  were  also  analyzed  using 
amplitude  data  alone „ 

Deisler  and  Wilhelm  (3)  also  employed  the  frequency  response 
technique  to  investigate  both  pore  diffusion  within  particles  and 
longitudinal  dispersion  in  the  fluid  between  particles  in  the  lami- 
nar  flow  of  gas  mixtures „   A  diffusion-type  model  was  used,  and  the 
analytical  solution  was  derived  based  on  B,C„(III)„   The  axial  dis- 
persion was  evaluated  simultaneously  with  the  diffusivlty  in  solid 
and  the  adsorption  constant „  Applying  a  similar  model,  boundary 
conditions  and  experimental  technique,  McHenry  and  Wilhelm  (13) 
studied  the  axial  dispersion  of  binary  gas  flowing  through  beds  of 
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inert  spherical  particles »   The  results  were  evaluated  using  am- 
plitude data  alone o 

A  very  interesting  paper  was  published  by  Head,  Hougen  and 
Walsh  (7)  in  their  study  of  the  axial  dispersion  of  water  flowing 
through  a  vertical  tubec  A  diffusion-type  model  was  assumed  to- 
gether with  BoCo(III)o   Experimental  data  were  obtained  by  the  pulse 
wave  technique,,  and  the  results  were  analyzed  by  reducing  pulse  data 
to  the  frequency  response  form  by  means  of  a  Fourier  transformation,, 
The  reduced  data  were  plotted  in  a  Bode  diagram  to  evaluate  the 
mean  dispersion  coefficient  by  fitting  experimental  data  and  calcu- 
lated values o  A  similar  method  of  analyzing  experimental  data  was 
reported  by  Romano  and  Pigford  (17)  in  their  study  of  the  perform- 
ance of  a  fluidized  bed  by  the  frequency  response  technique 0   The 
data  obtained  were  reduced  to  pulse  form  by  a  Fourier  transformation, 
and  the  resulting  residence  time  distribution  curve  was  used  to 
study  the  mechanism  of  gas  mixing  in  a  fluidized  bed,  and  the  effects 
of  gas  velocity  arid  bed  height  on  the  distribution  function,, 

Otake  and  Kunugita  (15)  presented  a  general  discussion  of  mix- 
ing  characteristics  of  reactors  without  chemical  reaction  based  on 
•a  diffusion- type  model0   Analytical  solutions  corresponding  to  the 
different  sets  of  boundary  conditions  were  presented  in  their  paper 
for  various  types  of  transient  inputs,,   In  another  paper  (16),  they 
investigated  the  performance  of  irrigated  packed  towers  using  a 
diffusion-type  model  and  B*C.(II)0  A  step  signal  was  used,  and  the 
data  were  analyzed  under  the  assumption  that  when  the  degree  of 
mixing  is  very  small  their  analytical  solution,  which  is  in  the 
form  of  an  n-th  term  expansion,  can  be  simplified  to  only  the  first 
term0 
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In  the  previous  section,  it  was  pointed  out  that  many  other 
mathematical  models  have  been  proposed  to  predict  the  behavior  of 
continuous  flow  systems,,   One  of  the  most  widely  used  models  be- 
side the  diffusion- type  model  Is  the  diffusion-dead  space  modelo 
This  model  for  axial  dispersion  postulates  the  existence  of  a  sta- 
tionary fluid  region  In  the  voids  of  the  process  equipment  which 
can  exchange  mass  with  the  moving  fluld0   In  a  purely  theoretical 
Investigation,  Turner  (20)  presented  a  model  of  a  packed  bed  which 
can  be  classified  under  the  diffusion-dead  space  category.   He  In- 
troduced a  concept  hydrodynamlcally  equivalent  to  the  packed  bed, 
which  consists  of  a  number  of  Identical  channels  through  which  the 
fluid  Is  flowlngo  Each  of  these  channels  has  a  number  of  "pockets" 
of  stationary  fluid  which  are  arranged  along  the  main  channel 0   The 
fluid  within  the  "pockets"  Is  regarded  as  stationary 0   The  mass 
transfer  between  the  flowing  fluid  and  the  "pockets",  and  the  mass 
transfer  within  the  "pocket"  volume  Is  assumed  to  take  place  by 
molecular  diffusion,.   The  model  proposed  by  Turner  (20)  appears  to 
be  too  complex  to  be  of  practical  use0   In  an  attempt  to  simplify 
the  model,  Koump  (9)  presented  the  diffusion-dead  space  model  under 
the  assumptions  that  the  axial  dispersion  of  mass  takes  place  by  a 
combination  of  two  mechanisms:   a  diffusion  type  process  and  a  mass 
exchange  between  the  moving  fluid  and  the  fluid  within  the  stagnant 
void  volume o   The  experiment  was  carried  out  by  three  different 
techniques  s   frequency  response,  Impulse,  and  step  Inputs  <>  The 
results  were  analyzed  by  both  the  diffusion -type  model  and  the 
diffusion-dead  space  model. 


13 


Von  Rosenberg  (23)  used  the  diffusion-dead  space  model  to  In- 
vestigate the  problem  of  the  displacement  of  a  fluid  from  a  bed  of 
inert  solid  by  another  fluid 0   He  regarded  the  packed  bed  as  an 
assembly  of  capillaries  and  used  the  solution  of  G.  I.  Taylor  (19) 
for  the  case  of  negligible  forward  and  rapid  lateral  diffusion  as 
a  basis  of  the  analysis 0   He  assumed  that  the  axial  dispersion  takes 
place  by  the  combined  effects  of  velocity  distribution  in  the  pores 
and  rapid  diffusion  in  the  radial  direction,.   The  concept  of  an  ax- 
ial dispersion  coefficient,  therefore s  does  not  enter  the  picture 
in  this  approach. 

Other  models  which  are  often  used  in  the  study  of  axial  mixing 
are  the  cell-model  and  cell  with  back-mixing0   The  concept  of  the 
cell-model  implies  that  the  spread  of  injected  tracer  as  a  result 
of  passage  through  a  bed  can  be  described  by  the  equation  of  con- 
centration-time behavior  in  the  effluent  of  the  n-th  perfect  mixing 
cello   The  cell  with  back-mixing  differs  from  the  cell-model  in  that 
each  one  of  the  n  perfect  mixing  cells  is  associated  with  a  definite 
back-mixing. 

Carberry  and  Bretton  (1)  investigated  the  longitudinal  disper- 
sion during  the  flow  of  water  through  beds  of  rings  by  the  cell- 
model.   Two  different  types  of  input  signals  were  used  in  their  work; 
an  impulse  function  which  injected  the  tracer  instantaneously,  and  a 
rectangular  input  function  which  injected  the  tracer  within  a  small 
given  time.   In  both  cases,  the  results  were  analyzed  by  fitting  the 
shape  of  the  experimental  response  curves  to  the  calculated  response 
curves  based  on  the  cell-model.   Miyauchi  (Ik)    presented  a  purely 
theoretical  analysis  of  the  cell  with  back-mixing  model. 
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It  seems  evident  from  this  literature  survey  that  the  diffusion- 
type  model  has  been  employed  almost  exclusively  because,  to  a  first 
approximation,  it  gives  a  reasonable  description  of  the  phenomenon 
associated  with  mixing  inside  reactors.  There  seems  to  be,  however, 
some  disagreement  on  the  method  of  analyzing  experimental  data  and 
on  which  sets  of  boundary  conditions  to  apply,,  All  of  the  mathemat- 
ical solutions  are  most  certainly  oversimplified  and  unrealistic 
from  the  standpoint  of  representing  the  true  mechanismc   Due  to  the 
simplification  involved  in  the  mathematical  treatment  on  B#C.(III), 
it  has  been  used  most  often  without  much  quantitative  justif ication0 
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THEORY 

The  purpose  of  this  section  Is  to  derive  the  partial  differ- 
ential  equation  for  the  diffusion -type  model  of  axial  dispersion 
considered  in  this  investigation,  and  to  present  the  mathematical 
expressions  for  four  sets  of  boundary  conditions c 

Diffusion  Model 

This  model  assumes  that  the  rate  of  dispersion  is  proportional 
to  the  bulk  concentration  gradient  of  the  diffusing  component  in 
the  tubular  reactor „  Under  this  assumption,  a  material  balance  with 
respect  to  a  traceable  component  in  the  fluid  with  chemical  reaction 
occurring  gives  the  following  conservation  of  mass  equation  in  vec- 
tor forms 

-|^  +  div(-D  grad0  C)  +  div(uC)  +  r(C)  =0  (1) 

where  D  is  a  dispersion  coefficient  or  apparent  diffusivity  which 
characterizes  the  dispersion  of  fluids  flowing  with  a  velocity  of 
u,  and  C  and  t  represent  the  concentration  and  time  variables  re° 
spectivelyc  The  first  term  of  equation  (1)  represents  the  change 
in  concentration  with  respect  to  time0  The  second  term  represents 
the  net  out-flow  of  component  due  to  dispersion,,  The  third  term 
represents  the  velocity  gradient  while  the  fourth  term  represents 
the  rate  of  chemical  reaction,, 

The  following  assumptions  are  made; 

1)  Only  the  dispersion  in  the  axial  direction,  z,  Is  signifi- 
cant,  that  is,  the  physically  three-dimensional  problem  of  the 
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second  term  in  equation  (1)  is  idealized  to  a  one-dimensional  trans- 
port equation,, 

2)  The  coefficient  D,  in  equation  (1)  is  independent  of  position 
and  concentration  gradients  as  in  the  case  of  fully  developed  turbu- 
lent flow,  and  is  characterized  by  D„ 

3)  Changes  in  volume  are  assumed  not  to  take  place,  so  that  the 
velocity,  u,  in  equation  (1)  is  the  same  at  all  cross-sections,  and 
is  characterized  by  u„ 

k.)    Only  the  chemical  reaction  of  first-order  is  considered,. 
With  these  simplifications,  equation  (1)  becomes  the  following  par- 
tial differential  equation? 

JL£  m  £  _d_c  _  -  £c   _  kc  (2) 

where  k  is  the  reaction  rate  constant,, 

Equation  (2)  is  a  second-order  linear  partial  differential  equa- 
tion which  describes  the  dispersion  phenomenon  in  terms  of  time,  po- 
sition and  dispersion  coefficient,,  As    stated  in  the  previous  section, 
the  solution  of  equation  (2)  depends  on  the  initial  and  boundary  con- 

ditions<> 

Boundary  Conditions 

Various  boundary  conditions  have  been  suggested  for  use  with 
equation  (2)0  As  stated  in  the  previous  section,  boundary  conditions 
may  have  a  significant  effect  on  the  solution  of  equation  (2)<> 
Hence,  the  boundary  conditions  used  with  equation  (2)  must  be  realis- 
tic and  quantitatively  justified.   One  of  the  purposes  of  this  inves- 
tigation is  to  find  boundary  conditions  compatible  with  the  solution 
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of  equation  (2) „ 

Danckwerts  (2)  introduced  a  set  of  intuitive  boundary  condi- 
tions based  on  the  equation  of  continuity 0   He  proposed  that,  owing 
to  the  dispersion  taking  place ,  the  concentration  of  reactant  in 
the  entering  stream,  C]_,  must  be  larger  than  the  concentration  just 
within  the  entrance  of  the  reactor. 
This  was  expressed  mathematically  as; 

U<J1=UCZ_0+  -B0!)z_o+  (3) 

At   the   outlet   of  a  finite   tubular  reactor,    he   proposed   the 
same   relationship  as    equation   (3)t 

uCz_L-  -  uCz_L+   -  d(-||)z_l+  (3a) 

The  equation  of  continuity;,  however,  necessitates  equation  (3a)  to 
be  expressed  ass 

uCz_L-  -  SCZ_L+  (3b) 

Substitution  of  equation  (3b)  into  equation  (3a)  yields; 


3z'z— L' 


)    T+  =  °  (3o) 


Since  D  is  different  from  zero,  the  term  in  parenthesis  of  equation 
(3c)  has  to  satisfy  the  following  expression; 
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Equations  (3)  and  (lj.)  will  be  designated  as  B.C. (I). 

In  an  attempt  to  approximate  the  above  B.C. (I),  Hulburt  (8) 
proposed  a  set  of  boundary  conditions  which  neglect  the  effect  of 
diffusion  at  the  reactor  entrance,,  Under  this  assumption,  the 
second  term  in  the  right-hand  side  of  equation  (3)  vanishes ,  and 
the  following  boundary  condition  at  the  inlet  of  the  reactor  can 
be  obtained? 

uCx  tt   uCz_^0+  (S) 

At  the  outlet,  z  =  L,  he  used  the  same  boundary  condition  proposed 
by  Danckwerts.   Equations  (5)  and  (lj.)  will  be  designated  as  B.C. (IT). 

Levenspiel  and  Smith  (12)  and  several  other  investigators  (ll|) 
(15)  proposed  a  set  of  approximate ■ boundary  conditions.   They  not 
only  neglected  the  effect  of  diffusion  directly  within  the  entrance 
of  the  reactor,  but  also  treated  a  finite  length  of  reactor  as  an 
infinite  reactor,  i.e.,  one  value  of  the  dispersion  coefficient 
applied  to  the  entire  length  and  cross-section.   This  infinite 
length  reactor  boundary  condition  is  expressed  mathematically  as; 

lim  C(z,t  )=  0  (6) 

Equations  (5)  and  (6)  will  be  designated  as  B.C. (III). 

Another  set  of  boundary  conditions  which  has  never  been  employ- 
ed previously  is  presented  in  this  investigation.   This  set  of  bound- 
ary conditions  preserves  continuity  just  as  those  of  Danckwerts . 
The  only  approximation  made  in  this  set  of  boundary  conditions  is 
that  the  finite  length  of  reactor  is  treated  as  a  section  of  an  in- 
finitely long  reactor  as  in  the  case  of  B.C. (III).  The  feasibility 
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of  this  set  of  boundary  conditions  will  be  discussed  in  a  later  see- 
tion  in  connection  with  steady-state  first=order  conversion  and  the 
response  to  various  transient  inputs 0   Equations  (3)  and  (6)  will  be 
designated  as  B.C. (IV), 

Table  2  summarizes  the  proposed  boundary  conditions  for  use  with 
equation  (2)0 

Table   2<>      Summary   of  proposed  boundary   conditions 0 


Boundary  Conditi 

ons 

Inlet    (z 

=  0) 

Outlet    (z  ~  L) 

B.C. (I) 

uCi 

=  uCz^o+ 

B.C. (II) 

uC;l 

=  ut:z^o+ 

(J£\             =  o 
vaz^z  .  L 

B.C. (Ill) 

UCX 

■  uCz_^o+ 

lim(z,t)    =  0 

Z.-+00 

B.C. (IV) 

uCi 

=  uCz— -G+ 

lim(zst)    =  0 

Z-+CO 
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MATHEMATICAL  TREATMEM1 

The  solutions  of  equation  (2)  which  satisfy  the  four  sets  of 
boundary  conditions  tabulated  In  Table  2  of  the  previous  section 
are  derived  by  means  of  Laplace  trans format Ions 0   This  has  the  ad- 
vantage of  transforming  a  partial  differential  equation  into  an 
ordinary  linear  differential  equation0   To  derive  the  expressions 
for  the  exit  concentrations  against  various  types  of  transient  in- 
puts,, the  back  transformation  Is  carried  out  for  each  of  the  solu- 
tions  by  means  of  tables  of  Integral  transforms  (5)  or  the  Heaviside 
partial  fraction  theorem,. 

This  section  will  be  divided  Into  four  parts  to  show  how  the 
solutions  are  obtained  for  both  ideal  and  non-ideal  flow  In  steady- 
state  and  unsteady-state  isothermal  tubular  flow  reactors 0   Only 
the  final  results  will  be  given  herep  and  the  details  of  the  deri- 
vations will  be  given  in  Appendix  I,  II,  III,  and  IV<, 

Steady-state  Non-ideal  Flow 

The  differential  equation  for  the  steady-state  tubular  flow 
reactor  based  on  the  diffusion- type  model  with  axial  dispersion  and 
first-order  chemical  reaction  can  be  obtained  from  equation  (2)  by 
letting  the  accumulation  be  zero? 

o  =  50  -  Saf  -  kc  (7) 

The  following  transformations  of  the  variables  in  Equation  (7) 
lead  to  the  dimensionless  form  of  the  differential  equation? 
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<*£.   -   2M  i^-   °   2MR  /    =0 
a  )72  a  >? 

(8) 
where         j/=  C  »  =  z    .    M  B  uL  and  R  =  kL 

cl  L  2D  u 

The  four  sets  of  boundary  conditions  tabulated  in  Table  2  of 
the  previous  section  are  transformed  likewise  and  are  summarized  in 
Table  3. 

Table  3.   Transformed  forms  of  boundary  conditions,, 


Boundary  Conditions     Inlet  (  ^  =  0)         Outlet  (  >[   =  1) 


B.C. (I)  yL   a  2M  (Y(0+)  -  1  ]        %L   =  0 

B.C. (II)  /  =  1  d/  _  0 

d^ 

B.C. (Ill)  ^  =  1  lim  /  =  0 

B°C'(IV)  |£  -  2M  [/(0+)  -  1  ]        W  =  0 


The    details    of   the   solutions    of  equation   (2)    which  satisfy   the 
above  boundary  conditions    (Table  3)    will  be  given   in  Appendix  Ic 
The   results    obtained   for  each   of  the   above   four  boundary  conditions 
are    respectively? 

yT(>7  s    _   2f(l+3)exp[M(l+<6)    +   M(l-0)fr]    =    (l-jfl)exp[M(l+^  +   M(l-6fll 
^l("]        2[  d+/3)aexp[MU+p)3     -    (1^/expcMU^)] ^J 

y>      (>9i    .    (l+l)exp[M(l+0)    +   M(l-0)£1      -    ( 1-/8)  exp[M(  1+0)  ?t    +   M(l-0fl 
^(  Y  } (l+p)expmi+p)l     -    (l-^expcMd-^)] ~ ,  ., 

»  — 


(9) 


0) 
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y__T(^  )  =  exp\M{l-p))Q  (11) 


III 

yiv(>?)  a  2exp[M(l^)>eJ  (12) 

where    /?  =  /l  +  *R 

Steady-state  Ideal  Plow 

In  practice  neither  plug  flow  nor  completely  mixed  flow  will 
occur;  there  always  will  be  some  degree  of  dispersion  due  to  viscous 
effects  and  molecular  or  eddy  dlffusione   In  order  to  establish  the 
theoretical  limitations  involved  with  regard  to  the  degree  of  mixing, 
It  Is  necessary  to  derive  the  equations  for  the  Ideal  flow  reactor0 

For  the  completely  mixed  flow  reactor,  where  D  =co   and  M  =  0, 
and  If  the  effective  volume  of  the  reacting  mixture,  V,  and  the  volu- 
metric flow  rate,  v,  are  constant,  the  material  balance  Is  the  same 
as  for  a  single  uniformly  stirred  tank  reactor  and  may  be  written  as: 

V  ^  =  vC1  -  vC  -  kCn  (13) 

where  n  Is  the  order  of  reaction,. 

For  the  calculation  of  the  steady-state  concentration  profile 
and  the  conversion  for  the  case  of  a  first-order  chemical  reaction, 
equation  (13)  may  be  simplifledo 

0  =  JL  Ci  -  -2-  C  -  £  C  (Ik) 

Equation  (11+)  can  be  transformed  Into  a  dlmensionless  form  by 
the  substitutions  used  for  (8)  to  obtain: 
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r/+/„1  =  0  (15) 


or 


For  the   steady-state  plug  flow  reactor,   where  D  -  0   and  M  -co  , 
the   first   term  in  the   left-hand  side    of  equation   (8)    vanishes    to  give 

f^-^-R^  (17) 

The  solution  of  equation  (1?)  is; 
l/iy  )   =  e-R*  (18) 

Because  of  the  theoretical  limits,  which  were  pointed  out  early 
in  this  section,  all  the  concentrations  and  conversion  calculated 
from  equations  (9)  through  (12)  have  to  lie  between  the  results  of 
equations  (16)  and  (18),  which  give  the  minimum  and  maximum  values 
respectively* 

Unsteady-state  Non-ideal  Flow 

When  the  deviation  of  the  concentration  of  the  traceable  com- 
ponent from  the  steady-state  value  due  to  a  pertubation  of  Inlet 
concentration,  C^(t),  only  is  considered,  the  initial  condition  may 
be  set  as  C(z,0)  ■  0o  The  transfer  functions,  W(s),  defined  as  the 
ratio  between  the  Laplace  transforms  of  the  inlet  and  the  outlet 
concentrations,  using  equation  (2)  for  each  set  of  four  boundary  con- 
ditions are  shown  in  equations  (19)  through  (22) „   The  details  of 
the  derivations  are  shown  in  Appendix  II 0 


2h 


W(s)  -  d+^expC-Md-^')]  -/li-^)*eXpc-Ml+^)]        U9) 

2.1s' 
W(3)  "  d+^)  expc-M(l-^)j  -  (Up)    expc-Md*^)]        (20) 

W(s)  a  expl>d-/S')]  (21) 

W(a)  =  ^xp[M(l-^)]  (22) 

r 

The  basic  forms  of  the  inlet  concentration  variations s  C^t), 
and  their  Laplace  transforms ,  C-j_(s),  for  impulse  and  step  input  ares 

C^t)  C-^s) 

Impulse  C^u'Ct)  C1 

Step  Input  C-ju(  t)  G^/a 


The  transforms  of  the  exit  concentration,  C2(L,s),  are  then  obtained 
systematically  as  the  products  of  the  transfer  functions,  W(s),  and 
the  C]_(s)  as  follows; 

Impulse  C2(s)   =   C]W(s) 

Step  Input  C2(s)   =   C-]W(s)/s 

The  inverse  transforms  of  C2(s)  to  obtain  the  response  of  exit 
concentration,  C2(t),  against  the  variation  of  feed  concentration, 
C^(t),  may  be  carried  out  systematically  by  means  of  the  Heaviside 
partial  fraction  theorem,,   The  results  are  tabulated  in  Table  l\.   for 
each  of  the  boundary  conditions 0   The  details  of  the  derivations  are 
given  in  Appendix  III* 
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Table  l\.„      Response    of  exlc   concentration  from  an    isothermal 
tubular   flow   reactor,, 


Transient 
Inputs 


B.C„(I) 


Impulse 


c2v 


oo 


2    Sn    (M  sin    Sn    +     Sn  cos   <fn  ) 

(M2  4-   2M  +    S^-    ) 


0s| 


x  Exp 


m  -  fyi2  ±  Sn2  *  am'i 

2M  ^ 


(23) 


c2(t) 

-  A    -  k 

CO 

c. 

1 

n= 

Step 

Input 

— 

x  Exp 

M 

M  J"^  (M  sin   ^    +    Jn  cos^fc) 
(M2  +   2M  +  <f  2)  (M2  +  <f  2  +   2MR) 


M  YM2  +^k 2  +   2MR\ 


•Q 


*2lT 


J€ 


(27) 


Frequency 
Response 


A*    R.    = 


4>    - 


b   e 


M 


z 


(  «<x  +    ^<  2) 


m      -1/    <K  2x 
Tan  x( — li)    +   a 

<*1 


(3D 


Table  2j.„   (Cont0) 
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Transient 
Inputs 


B0Co(II) 


Impulse 


c2v 


oo 


-  sn 


-^7      (M+   M2  +     ^2}    cos  ^' 


x  Exp      M  «/M2  +  <£>'2  +   2MR  \  /r 

I 2w— — -v  6 


(2i|) 


Step 
Input 


Q? 


=     D     •"     C 


1 


^77  (M2  +  S^2*  M)  (M2  +^2+2MR) 


x  Exp 


M 


M2  + 


<£/2  +   2W 

"2H 


> 


(28) 


Frequency- 
Response 


M 


Ao    R, 


b   e 


1     '  M1  ■)      ■[ 

(    ^3  +     oik)* 


4>      =  a   -  Tan"1   ( °fr-      ) 

<*3 


(32) 


Table  k0   (Cont.) 
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Transient 
Inputs 


B.C. (Ill) 


Impulse 


c2v 


.-1.5 


x  Exp 


M   -  Jr(    M£  +  £  +    2RC  ) 


(25) 


c2.(t) 
"K 


_  i 

—     2 


M+(M2  +    2MR)2 


Step 
Input 


x  Erfc[|^)S-J-|(2M€+  ljBO* 
+   eM"(M2  +   2MR)2 


x  Erf  c 


*6r")-  ir(2M^  +  UE€)* 


€ 


Frequency 
Response 


Ao  Ro  =  eM"x 


4)    =- 


y 


(33) 


Table  lj.0   (Cont«) 
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Transient 
Inputs 


B.C. (TV) 


Impulse 


C2V   _/2M^2  f/r\  "*  x  Exp|M-|(MC    +  H  +   2R  € ) 


-M  Exp(2M  -  R^)Erfc|  f  Jl  \ 


(26) 


Step 
Input 


c2(t) 


Jo 


(2M](X)"*     Exp   M-|(M\+  |  -   2R\  )' 


-M  Exp(2M  -  Rt)    Erfc 


fef  *  (B* 


dA 


(33) 


Frequency 
Response 


A.    R,   = 


2   e 


M-x 


[(1  +   x/M)2+  sin2  a]  ^~ 


<p       =  -y   -  Tan-l^aln  a    \ 

•:-  x/M/ 


(314-) 
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The  amplitude  ratios  and  phase  shifts  of  the  frequency  response 
functions  caused  by  a  sinusoidal  Input  of  amplitude  C-j_  and  angular 
frequency  60  ,  are  obtained  by  means  of  a  Fourier  transformation 
Thus,  a  performance  function,  defined  as  the  ratio  of  the  Fourier 
transform  of  the  output  to  that  of  the  Input,  Is  obtained  by  replac° 
Ing  s  by  j£0  In  equations  (19)  through  (22)  for  each  of  the  boundary 
conditions  respectively,  where  j  =  /^T.   From  these  performance  func< 
tlons,  the  computable  forms  of  the  amplitude  ratios  and  phase  shifts 
are  obtained  by  changing  Into  polar  coordinates »   These  are  also 
summarized  In  Table  fy.«   The  details  of  derivations  will  be  given  In 
Appendix  IT, 

Unsteady-state  Ideal  Flow 

For  a  completely  mixed  reactor,  where  D  =  oo     and  M  ~  0,  the 
Laplace  transform  of  equation  (13)  with  respect  to  time  with  zero 
initial  concentration  gives  a  transfer  function  of  the  following 
form: 

W  (s)  "  1  +  ia+k)e-  (35) 

Substitution  of  the  corresponding  transient  Inputs  and  their 
Laplace  transforms  into  equation  (35)  and  inverse  transformation 
of  the  resulting  equations  gives: 

Po.'t-Ur    -U+R)*7 
Impulse  igj^Jv  „  e  (35) 

r  -d+iO'C- 

(37) 


Step  input         £gl£i  =  ^  [l-e       J 
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Frequency  Response 


(1+R)2  +  ( t>)Q  )2 


Cp-f3m=o  (jc0)  = „ x     -  ;  e  (38) 


where   P8F0   ~  the  performance  function,, 

Q  =  quantity  of  tracer  injected  in  the  impulse  functionc 
V  ss  effective  volume  of  reactor 

/C=  t/e 

For  the  unsteady-state  plug  flow  reactor,  where  T)  =  0  and 
M  -co,    the  first  term  on  the  right-hand  side  of  equation  (2) 
vanishes,  so  that 

£2-   =  -u  -<L2.  ~  kC  (39) 

at       dz 

The  Laplace  transform  of  equation  (39)  with  respect  to 
time  with  zero  initial  concentration  gives  a  transfer  function 
of  the  following  form? 

-R-s© 
W(s)  «  e  (24-0) 

Substitution  of  the  corresponding  transient  inputs  and 
their  Laplace  transforms  into  equation  (2+0 )  and  inverse  trans- 
formation of  the  resulting  equations  gives t 

C?  (t)  V     -R    c  ~ 

Impulse  -S -a      d    (  C  -  €0  )  (kD 

Q 
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Step  Input 


C2    (t) 


0  0  <C<      1 

e"R  Z  >      1 


(1*2) 


_         __  -R     -j(A>0) 

Frequency  Response  ]_P8F0J  M  (jA))    =  e        e  (1^3) 
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NUMERICAL  COMPUTATION 

The  numerical  computations  for  steady-state  solutions  repre- 
sented by  equations  (9)  through  (12),  (16)  and  (18)  are  straight 
forward  and  those  for  the  frequency  response  solutions  represented 
by  equations  (31)  through  (3h)  >    (38)  and  (lj.3)  are  easy  to  perform 
by  use  of  a  computer,  and  no  elaboration  is  required.   This  section 
is  primarily  concerned  with  the  computations  for  unsteady-state 
transient  inputs,  impulse  and  step  inputs,  represented  by  equations 
(23),  (21;),  (27)  and  (28).  All  of  these  four  equations  are  in 
series  form  which  requires  that  a  tolerable  error  be  established 
such  that  the  error  term  of  the  series  does  not  exceed  10~-^« 

The  first  step  in  calculating  numerical  values  of  these  equa- 
tions was  to  determine  the  method  of  convergence  of  the  equations; 


cot 


S  =  Kf  "  t)       tor    ^  m 


tan  S     =  +fl-  for    $'n  (1+5) 

For  this   purpose,    the   Newton-Raphson   iterative   formula  was   used, 
name ly , 


tf"3  -  S*  -  x^i 


sc$J 


(l|-6) 


=     £~CfJ  +        -ct    SnP  -   k        K  ^ 


Jn 


6n    '  n  m      sfli*) 
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with  a  first  approximation  of  the  following: 


cCoJ 


■Co] 


Co] 


C-^J  cLOj 

6n     =    dn-i 


(kl) 


rr 


n  =  2,   3,   1+., 


•  •  •  e   « 


Equation  (lj.7)  holds  for  all  values  of  M, 

The  n-th  term  of  equations  (23),  (2lj.),  (27)  and  (28)  were  then 

calculated  with  the  d     and  S    thus  obtained  from  equations  (kl±)    and 

([(.5)  j  and  the  infinite  series  is  truncated  after  the  n-th  term  such 
CoV  x   /  n 


that 

error  of  10 


i 


n 

-5 


n=i 


c2v 

0. 


n 


is  smaller  than  the  maximum  tolerable 


The  block  diagrams  and  computer  programs  for  use  with  the  IBM 
65>0  computers  are  presented  in  Table  5>  for  equations  (23),  (2l|_)  and 
(1+4)  and  in  Table  6  for  equations  (27),  (28)  and  (1+5)  for  R  =  0 
only. 
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Table    5  -a       Block    diagram     for     use     with     I B  M    650 
for     computing       equations     (23),    (24)     and      (44). 


START 


R  CD(M)in  1951 


I 


Set    I.  R.A.=-8 


Trans  1959  A  to  1977    and 
check   M^O    if  M=0— *-halt 


F*- 


j 


Set     -t  =  <x 

*X.  increment  —  Vt  increji 


CQ3 


Setin=t^2,     n=l,    J^-=o 


Calc.in       ,    equ.   (4  6) 


fo+il 


<£ 


No 


-sr 


/tq] 

£n  ' 


ztq+0 

\?n 


^Yes~ 


Calc.   equ.  (23)  or  (27) 


€. 


yl/Yes 
Punch    M,  N,  <,~<r 


Add    I   to   I.R.A 

~ — IT     " 

.(l.R.A  =0  ) 
Yes 


<^-^+< 


increment 


Resit.   R, 


35 


Table  5-b,  Computer  programs  for  use  with 
equations  (23)  (24)  and  (44) 


IBM  550  for  computing 


MOD  IF 


LDD 

STD 

32 

LDD 
STD 


1582 
1570 
1500 
1563 
1571 


0700 


DEL17 


1563  -  Number   of  R 

Program  Step  No.  202,  Connect  to  1570 
FAD  1500  0700 
RAU  1570 
A  UP  1590 
STU  1570 
RAU  1571 
SUP  1590 
STU    1571 

NZU    NO       BACK  to  prog, 
1590  =  Fixed  point  one 
R 
R 
R 
R 
R 
R 
RSA    CCO  8     MODIF 
Program  step  No.  255,  Connect  to  1560 


1572 
1502 
1532 
1553 
1573 
1563 
0277 
1570 
1560 
1580 
1551 
1561 
1581 
1552 
1552 
1590 
1500 
1501 
1502 
1503 
1504 
1505 
0151 
0^27 


59  1582 
24  1570 
32  1500 
69  1563 
24  1571 
00  0005 
34  0250 
32  1500 

60  1570 

10  1590 
21  1570 
60  1571 

11  1590 
21  1571 
44  0030 
CO  0001 
00  0000 
10  0000 
10  0000 
50  0000 
3.0  0000 
50  0000 
81  0008 
71  1877 


1592 
1553 
0700 
1573 
0107 
0000 
1570 
0700 
1580 
1551 
1561 
1581 
1552 
1552 
0777 
0000 
0000 
0050 
0051 
0051 
0052 
0052 
1572 
1560 
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r  if  pot 

8  TO 

E  X  P03 

15 

0000 

2  4 

00  0  3 

0  006 

8  TU 

EXP03 

16 

000  6 

_._?i 

00  10 
80  0  7 

001  J 

LOO 

8  007 

17 

0  0  13 

69 

0019 

S  TO 
R  A  C 

E  X  P04 

18 
19 

0019 
00  3  5 

24 
68 

00  2  2 
000  0 

002  5 

00  00 

003  1 

R  A  U 

E  X\P03 
E  X  P05 

EXP  06 

•  30 
2  1 

0031 
0015 

60 
46 

0010 
00  18 

00  15 

BM  1 

006  9 

-LXJ--QjS_ 

F  BR 

LN.T..EJ. 

__2_2-_ 

-  0069 

33 

46 

00  7  2 
00  0  3 

00  4  9 

BM  1 

E  X  P07 

23 

00  4  9 

005  3. 

_   A  XC 

000  1 

E  XP  06 

24 

005  3 

58 

000  1 

006  9 

E  XP07 

F  A  0 

L  N  T  EN 

EXPOS 

35 

000  2 

32 

00  7  2 

009  9 

EXPOS 

3  X  C 
F  A  0 

0001 

36 

27 

00  18 
00  2  4 

59 

33 

000  1 

00  7  3 

0  03  4 

LN  TEN 

014  9  ' 

BM  1 

EX  EQ5 

EXP  08 

28 

0  14  9 

46 
21 

0018 
00  10 

009  9 

EXPOS 

8  TU 

E  X  P03 

29 

0099 

006  3 

S  T  U 

E  X  P09 

30 

00  6  3 

00  2  1 

31 
65 

00  6  8 
00  2  2 

00  2  1 

R  AL 

E  X  P04 

31 

002  7 

i    loo 

8  TO 

8  007 
E  X  P04 

.  32 

33 

00  2  7 
003  3 

69 

2  4 

80  0  7 
00  2  2 

0033 

007  5 

R  A  C 

8  0-0.3 

34 

007  5 

88 

80  0  2 
00  3  6 

0083 

LOO 

FPONE 

35 

008  3 

69 

00  3  9 

8  TO 

E  X  P10 
E  X  PI  1 

E  X  P  1  2 

36 
37 

0039 
00  4  5 

24 
34 

00  4  2 
00  4  8 

00  4  S 

S  TO 

000  1 

EXP18 

..  R  A  U 

'  E  X  P10 

38 

000  1 

60 

00  4  2 

004  7 

FA  0 

EXPO? 

39 

004  7 

33 

006  8 

009  5 

Jill. 

EX  P-lp 

_  4D. 

009  5 

21 
60 

00  4  2 
00  4  8 

014  5 

R  A  U 

E  X  Pll 

41 

0  14  5 

010  3 

F  A  0 

FPONE 

4  3 
4  3 

010  3 
0113 

32 
21 

00  3  6 
00  4  8 

0113 

8  TU 

E  X  Pll 

005  1 

R  A  U 

E  X  P  0  9 

'   

4  4 
4  5 

0051 
0023 

60 
39 

00  6  8 
00  10 

002  3 

FMP 

EX  P03 

00  60 

^-0  V 

-  E-  X  P  11.  - 

4  6 

00  6  0 

3  4 
31 

00  4  8 
00  6  8 

00  9  8 

STU 

E  X  P09 

47 

009  6 

007  1 

F  0  V- 

EX  P 10  - 

.48 

00  7  1 
00  9  2 

34 
3  3 

00  4  2 
019  5 

009? 
012  1 

F  SB 

SIZES 

4  9 

BM  1 

EiP12 

-SO 

0  121 
007  4 

4  6 
60 

00  7  4 
0042 

000  1 

R  A  U 

E  X  P  1  0 

51 

009  7 

_ LUP 

10 

E  X  P04  - 

EXPO? 

52 

009  7 
00  36 

10 

10 

002  2 
00  0  0 

00  0  3 

~f"p~one~ 

0000 

0051 

53 

0  0  5  1 

-SIU.&- 

LO- 

0  0  0  0  . 

"04  1 

—  5  4 

019  5 

10 

00  0  0 
02  5  8 

00  4  3 

LNTEN 

23 

02  58 

5  1  SI 

55 

007  3 

33 

515  1 

F  nn  c  B 

JLLD.  E_l 

n  fc  i<n 

C_0_3 

5-6 

0050 
00  5  6 

24 
4  5 

015  3 
0110 

00  5  6 

NZE 

C  S  V03 

C8  V  04 

ooa 

57 

0011 

C  S  V04 

R  A  U 

FPONE 

C  3  V  03 

003 

---58 

00  11 

60 

00  3  6 

015  3 

C  8VO  3 

BM  1 

C  S  V  OS 

C  8  V  06 

004 

S9 

0110 

46 

016  3 

00  14 

C  8  VO  5 

-E-A-D- 

T  IT  0  P  1 

£05 

6  0 

0163 
0  0  4  3 

32 
4  5 

00  16 
00  4  6 

0  0  4  3 

NZE 

CSV  04 

006 

61 

0011 

RM  1 

C  S  V  0  5 

007 

62 

00  4  6 
0100 

46 
33 

016  3 
030  3 

0100 

FSB 

ON  E  P  1 

008 

63 

00  2  9 

BZE 

0  a  V  O7 

c«»nn 

0  0  9- 

64 

00  29 

45 
61 

00  3  2 
00  3  6 

0  13  3 

CS  V08 

R  S  11 

FPONE 

C  3  V  03 

0  10 

65 

0  13  3 

015  3 

C  S  V  0  6 

-FSB 
NZE 

T  *  0  P  1 

Oil 

66 
67 

0014 
009  3 

33 

45 

00  16 
00  9  6 

009  3 

CSV  04 

012 

00  11 

BM  1 

nuns   , 

0  13 

66 

00  9  6 
0199 

4  6 
32 

0  19  9 
020  3 

00  14 

FAD 

ON  E  P  1 

014 

69 

007  9 

»H 

fSini 

r.  a  v  nn 

(11  5 

_  _7-D„ 

.0079 

4  5 

00  3  3 
00  8  6 

013  3 

C  8  V07 

STU 

C  8  V09 

0  16 

71 

0033 

31 

00  8  9 

RSU 

FPONE 

0  17 

7  2 

0089 

61 

00  3  6 

00  4  1 

8TU 

C  S  V  10 

0  18 

73 

00  4  1 

21 

0  14  6 

024  9 

ST  U 

C  S.VL1 

0-19. 

7  4 

03  4  9 

21 
20 

00  0  4 
006  1 

000  7 

STL 

c  s  via 

CS  V  13 

0  20 

7  5 

000  7 

00  6  4 

fnnc) 

_SXO_ 

r  a  v  n  a 

A  F  A  1  N 

a  1  N 

7_6- 

0  1  5  0 

34 

0  15  3 
016  0 

010  6 

NZE 

C  8  V  03 

031 

77 

010  6 

45 

0153 

BM  1 

F  A  0 

C  S  V14 

C3  V  IS    

0  3  2 

7  B 
79 

0  16  0 
02  13 

46 
32 

0213 
00  16 

0  114 

CSV14 

T  t  OP  1 

0  2  3 

014  3 

NZE 

_JB  8  At  0  2 

0  3  4 

60 

014  3 
0196 

4  5 
46 

0  19  6 
0  2  13 

015  3 

8M  1 

C  S  V  14 

0  25 

8  1 

030  0 

F  SB 

ON  F  P  1 

0  3--L-  B_2 

0300 

53 

020  3 
00  8  2 

012  9 

NZE 

C  S  V16 

C  3  V  03 

0  27 
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a .     Block     diagram     for     use     with 
for     computing      equations     (27) 
and    C28) . 
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Table      6 
IBM    650 


b.     Block     diagram     for     use     with 
for     computing      equation     (45). 
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Table    6  — c       Computer     program      for      use      with 
IBM     650       for       computing .     equations      (27)    8  (45) 
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0  29  3 

69 

800  7 

00  4  9 

STO 

E  X  P04 

0  19 

207 

004  9 

24 

0  13  8 

0  14  1 

R  A  C 

8  002 

0  20 

208 

014  1 

88 

800  2 

009  9 

LOO 

F  P  ONE 

021 

20  9 

0099 

69 

004  0 

0  3  4  3 

S  TO 

E  X  P  1  0 
E  X  P  11 

E  XP12 

022 
0  23 

2  10 
211 

0  34  3 
014  9 

2  4 
24 

014  6 
010  2 

0  14  9 

STO 

0  3  55 

EXP  1Z 

R  A  U 

E  X  P10 

02  4 

212 

0  3  5  5 

60 

014  6 

0  30  1 

F  A  0 

E  X  P09 

0  25 

213 

0  20  1 

32 

0  3  3  4 

0111 

STU 

E  X  P10 

0  2  6 

214 

0111 

21 

014  6 

019  9 

R  A  U 

E  X  P  1  1 

0  27 

215 

0  199 

60 

010  2 

0057 

E«a. 

/PONE  _ 

028 
0  29 

216 
217 

00  5  7 
0  36  7 

32 
21 

00  4  0 
0  10  2 

0  3  6  7 

8  TU 

E  X  PI  1 

04  0  5 

R  A  II 

£  X  P09 

0  30 

21  B 

0  40  5 

60 

033  4 

013  9 

FMP 

E  X  PO  3 

031 

219 

0139 

39 

012  6 

0  17  6 

F  0  V 

E  X  P  11 

0  3  2 

220 

0176 

34 

010  2 

0  15  2 

STU 

E  X  P09 

03  3 

22  1 

0  15  2 

31 

03  3  4 

0  2  8  7 

fa.v 

E  X  P  10 

0  3  4 
0  35 

222 
22  3 

0  28  7 
0  19  6 

3  4 
3  3 

014  6 
00  7  6 

0  19  6 

FSB 

S  1  7.  E  8 

060  3 

BM  1 

E  XP  12 

0  3  6 

22  4 

0  60  3 

46 

0  15  6 

03  5  5 

R  A  U 

E  X  P  10 

0  37 

225 

0  15  6 

6  0 

0  14  6 

02  5  1 

AUP 

E  X  P04 

E  X  P  0  2 

0  38 

226 

0  2  5  1 

10 

0  13b 

Oil", 

F  P  ONE 

10 

0000 

OOS  1 

0  39 

227 

00  4  0 

10 

00  0  0 

00  5  1 

-    8  1  ZE8 

10 

-  0000 

004  3 

0  40 

228 

007  6 

10 

000  0 

004  3 

LNTEN 

23 

02  SO 

5  151 

0  41 

229 

0  188 

23 

02  5  8 

515  1 

Table      6-d       Computer        program       for      use      with 
IBM     650       for      computing       equations     (28)     8   (45). 


k2 


BLR 

1000 

19  99 

i 

OOOO 

00 

OOOO 

OOOO 

PR  E  S  1 

1  0 

OOOO 

00  4  8 

2 

OOOO 

10 

OOOO 

00  4  8 

1  N  C  D  T 

1  0 

00  0  0 

00  50 

3 

00  5  0 

10 

OOOO 

u  0  SO 

DE  C01 

1  0 

00  00 

00  50 

4 

0  10  0 

10 

OOOO 

00  5  0 

H»FPI 

15 

7  07  9 

6  3  5  1 

5 

0  15  0 

1  5 

70  7  9 

6  3  5  1 

F  P  0  4  3 

10 

Q  0  0  Q_  __ 

__  0  Q  4  3 

6 

0  20  0 
0  25  0 

10 
10 

OOOO 
OOOO 

00  4  3 

A  C  UC  V 

1  0 

00  00 

00  4  4 

7 

00  4  4 

1STAU 

5  0 

OOOO 

00  4  9 

8 

0  300 

50 

OOOO 

0  0  4  9 

M  X  T  ft  II 

36 

00  no 

00  5  1 

9 

0  J  5  0 

36 

OOOO 

00  5  1 

D   T  »  U 

5  0 

00  00 

004  9 

1  0 

0  4  00 

50 

OOOO 

00  4  9 

D  1  VGE 

9  9 

99  99 

99.9  9 

1  1 

0  4  50 

99 

9  9  9  9 

999  « 

ONE   5 

10 

5  0  00 

OQ  5  1 

1  2 

0  SOO 

1  0 

50  0  0 

00  5  1 

F  P   3 

30 

OOOO 

00  5  1 

1  3 

0  5  50 

30 

00  0  0 

00  5  1 

START 

R  C  0 

19  5  1 

1  4 

0  60  0 

70 

19  5  1 

000  1 

R  ft  U 

1  S  T  A  U 

SET 

1  5 

0  00  1 

60 

0  30  0 

000  5 

SET 

S  TU 

197  7 

1  6 

000  5 

21 

197  7 

00  30 

LOO 

0  E  C  01 

1  7 

00  30 

69 

010  0 

000  3 

.   S  TO 

1  N  C  0  T 

1  8 

000  3 

2  4 

00  5  0 

00  5  3 

L  0  D 

19  5  1 

1  9 

00  5  3 

6  9 

19  5  1 

000  4 

S  TO 

198  4 

2  0 

000  4 

24 

1  9  B  4 

00  37 

LOO 

F  P  0  N  F. 

2  1 

00  3  7 

69 

00  4  0 

004  3 

S  TO 

19  7  9 

A  K  1  0  » 

2  2 

00  4  3 

24 

197  9 

00  3  2 

ft  k  I  n  n 

R  A  II 

197  9 

2  3 

00  3  2 

60 

197  9 

00  3  3 

FMP 

ON  E  P  1 

2  4 

00  3  3 

39 

00  3  6 

0086 

STU 

N       P  1 

25 

008  ft 

21 

00  9  0 

00  9  3 

R  A  U 

197  9 

26 

009  3 

60 

197  9 

008  3 

FMP 

fPHO 

2  7 

008  3 

39 

0  13  6 

0186 

FSB 

F  P  ON  E 

28 

0186 

33 

00  4  0 

00  17 

FMP 

H  A  F  P  1 

29 

00  17 

39 

015  0 

06  50 

FAD 

1  N  C  0  T 

30 

0  65  0 

32 

00  5  0 

00  2  7 

FSB 

N   P  1 

31 

00  27 

3  3 

00  9  0 

HO  6  7 

N  2  U 

RE  T  U  N 

32 

00  6  7 

4  4 

00  2  1 

002  2 

BM  1 

RE  T  UN 

3  3 

00  2  1 

46 

00  2  4 

002  2 

FAD 

N   P  1 

34 

00  2  4 

5  2 

00  9  0 

0117 

STU 

L  0  C  T  0 

COttPU 

35 

0  117 

21 

00  7  2 

0  0  2  5 

C  OMP  u 

FD  V 

19  51 

36 

00  2  5 

34 

195  1 

00  5  1 

STU 

TEMPI 

37 

00  5  1 

31 

00  0  6 

000  9 

R  A  U 

LO  C  TO 

38 

000  9 

60 

00  7  2 

007  7 

LOO 

E  00  C  B 

39 

007  7 

69 

00  8  0 

0  13  3 

STU 

COS  L  N_ 

4  0 
4  1 

0080 
00  8  7 

31 

60 

00  3  4 
00  7  2 

00  97 

R  A  U 

L  0  C  TO 

0  12  7 

L  n  o 

E  00  C  A 

-  .-4  2 

0  12  7 
0  130 

69 
21 

0  13  0 
00  3  9 

018  3 
0137 

FD  V 

CO  S  1  N 

4  4 

0137 

3  4 

00  3  4 

0  134 

F  A  0 

TEMPI 

4  5 

0  13  4 

32 

000  6 

02  3  3 



BM  1 

ft  S  K  N  G 

CHECK 

4  6 

0  23  3 

46 

02  3  6 

0  18  7 

ISKIIG 

FAD 

P  R  E  S  1 

4  7 

0  2  3  6 

32 

OOOO 

0  17  7 

-  -  NiU 

N  0  3  1  R 

OK 

4  8 

0  17  7 

■  4  4 
46 

00  3  1 
018  4 

00  8  2 

NO  S  1  R 

RM  1 

MO  OF  V 

OK 

4  9 

0031 

OO  8  3 

MO  OF  Y 

R  A  U 

L  0  C  TO 

50 

0  18  4 

60 

00  7  2 

02  27 

F  A  0 

1  N  C  0  T 

51 

0  2  2  7 

32 

00  5  0 

027  7 

STU 

L  0  C  TO 

C  0  MP  U 

52 

027  7 

21 

007  2 

002  5 

CHECK 

FSB 

PR  E  SI 

5  3 

0  18  7 

33 

OOOO 

03  3  7 

HI  U 

T  E.  S  T  0 

n* 

5  4 

0  337 
008  1 

4  4 
46 

00  8  1 
00  8  2 

00  8  2 

TE  ST  0 

RM  1 

OK 

R  E  T  UN 

5  5 

00  3  2 

R  E  TU  N 

R  A  U 

1  N  C  0  T 

56 

00  2  3 

60 

00  5  0 

00  5  5 

F  0  V 

FP  T  »0 

57 

00  5  5 

34 

013  6 

0  2  8  6 

FSB 

F  P  04  3 

58 

0  28  6 

33 

02  0  0 

0  3  7  7 

RM  1 

MINUS 

59 

0  37  7 

46 

01  H  0 

013  1 

_   F  A  0 

F  P  0  4  3   . 

60 

0131 
0  4  27 

32 
21 

030  0 
00  S  0 

0  4  2  7 

STU 

1  N  C  0  T 

6  1 

010  3 

R  A  II 

L  0  C  TO 

63 

0  10  3 

60 

00  7  3 

0  4  77 

FSB 

1  N  C  0  T 

6  3 

0  47  7 

3  3 

00  5  0 

0  52  7 

STU 

LO  C  TO 

C  0  M  P  U 

64 

0  527 

21 

007  2 

002  5 

OK 

R  A  U 

L  0  C  TO 

65 

008  3 

60 

007  3 

057  7 

F_*LP_ 

LO  C  TO 

...68 

0  57  7 
0  12  2 

39 

21 

007  2 
00  2  6 

012  3 

STU 

SO  R   0 

67 

00  2  9 

R  A  U 

1951 

68 

00  29 
0  1  0  S 

60 
39 

19  5  1 
195  1 

010  5 

0  15  5 

STU 

M  SOR  E 

70 

015  5 

31 

00  10 

0013 

F  A  0 

S  0  R   0 

71 

00  13 

33 

00  2  6 

ol  5  3 

FJ1  V  . 

19.1L 

7? 

.015  3 

34 
34 

19  5  1 
013  6 

010  1 

FO  V 

FP  T  »0 

7  3 

0  10  1 

03  36 

FMP 

197  7 

74 

0  336 

39 

19  7  7 

06  2  7 

STU 

BO  X   3 

75 

0  627 

31 

013  2 

003  5 

R  A  U 

1951 

7  6 

00  3  5 

60 

19  5  1 

020  5 

FSB 

BOX   3 

77 

020  5 

33 

013  2 

O059 

LOD 

Eo_oj_a ._  . 

IB. 

00  59 

69 

00  12 

00  15 

STU 

~~ 198  2  ~ 

7  9 

00  12 

21 

19  8  2 

008  5 

9  A  U 

1951 

80 

008  5 

60 

19  5  1 

02  5  5 

FAD 

M  SOR  E 

81 

0  25  5 

32 

00  10 

02  37 

F  A  0 

SOR   0 

83 

03  37 

33 

00  2  6 

030  3 

STU 

BOX   3 

83 

0  30  3 

21 

000  8 

0011 

RAJ) 

JLtLC.IO- 

84 

0  0  11 

60 

007  2 
013  6 

0677 

FMP 

FFTIO 

85 

067  7 

39 

0386 

FMP 

SIGN 

86 

0  38  6 

39 

00  3  9 

008  9 

FO  V 

BOX   3 

87 

0089 

34 

000  8 

0058 

FMP 

1983 

88 

00  58 

39 

198  3 

018  2 

STU 

1981 

89 

018  2 

21 

198  1 

02  34 

FAD   SUM 

2  0 

0  33  4 
006  3 

32 

21 

03  8  7 
03  8  7 

006  3 

STU 

SUM 

91 

0140 

R  AU 

1981 

93 

014  0 

60 

198  1 

0135 

FOV 

SUM 

93 

013  5 

34 

02  8  7 

03  37 

RAM 

8003 

94 

0  337 

67 

800  3 

004  5 

R  A  U 

8003 

95 

00  4  5 

60 

80  0  2 

0253 

.  FSB 

A  C  II  C  Y 

?6_. 

0353 

33 

46 

02  5  0 
02  3  0 

07  2  7 

BM  1 

THANK 

97 

0  72  7 

018  1 

R  A  U 

197  9 

98 

0  181 

60 

19  7  9 

028  3 

F  10 

F  P  ON  E 

99 

0  38  3 

38 

00  4  0 

0167 

STU 

197  9 

100 

0167 

21 

197  9 

0333 

LOO 

OE  C01 

101 

03  3  2 

69 

010  0 

030  3 

SJ  Q_ 

1JLC  CL1 

.jki ao 

1.0  3_ 

.0303 

24 
69 

00  5  0 
00  7  2 

00  3  2 

TM  INK 

LDO 

LO  C  TO 

103 

02  30 

007  5 

STO 

1  9  8  0 

104 

00  7  5 

24 

198  0 

o  3  3  i 

43 


LOO 

S  U  GN 

105 

0  3  3  3 

S  T  0 

19B  3 

10  6 

0  387 

R  A  U 

V  P  0  N  E 

107 

0  43  6 

FSB 

SUM 

1  0  8 

0  0  9  5 

S  T  U 

197  8 

109 

0  113 

PC  H 

19  7  7 

A  0  T  A  U 

110 

0231 

»  n  t  «  u 

LOO 

19  99 

111 

077  7 

S  TO 

SUM 

11  2 

00  0  2 

R  A  U 

19  7  7 

113 

0  190 

FAD 

0   T  A  U 

. 

114 
115 

0  2  8  1 

FSB 

ONE   5 

0  8  2  7 

Bll  1 

GO  1 

116 

0  67  7 

F  A  U 

ONE   5 

SET 

117 

0280 

G  0  1 

FAD 

ONE   5 

118 

0  3  3  1 

FAD 

D   T  A  U 

119 

0  92  7 

_  FSB 

M  X  I  A  U 

130 
131 

0  97  7 

BM  1 

C  0  NN  1 

00  28 

•  F  AO 

M  X  T  AU 

SET 

123 

0  38  1 

C  ON  Nl 

LOO 

D  E  C  0  1 

133 

0  28  2 

S  TO 

1  N  C  0  T 
0  1  VG  E 

START 

124 
13  5 

0  3  5  3 

MINUS 

L  DO 

0  18  0 

STI). 

-1983 

_  P  UJM.C  K  

126 
127 

0  40  3 

PUN  C  H 

PC  H 

197  7 

0  4  6  6 

LOO 

19  99 

128 

0078 

S  TO 

19  8  3 

A  0  T  A  U 

129 

00  5  2 

E  OOC  B 

S  TO 

C  S  V  02 

BEC 

IN   COS 

130 

013  3 

NZE 

c  s  »n  3 

C  S  V  04 

002 

13  1 

0  139 

C  S VQ4 

RA.U 

_FP  ONI 

C  S  V  0  2 

0  03 

13  2 
1  3  3 

0  14  3 

C  S  V03 

BM  1 

C  S  V  OS 

C  S  V  06 

0  04 

00  4  2 

C  S  V  0  S 

FAD 

T  It  0  P  1 

0  05 

13  4 

014  5 

NZ  E 

C  S  V  04 

0  06 

135 

012  5 

BM  1 

c  s  vos 

007 

136 

0  128 

FSB 

0  N  E  P  1 

0  08 

137 

0  3  3  2 

. .  

NZE 

CSV  0-7 

C  S  V .0  8. 

009 

138 

0  16  3 

C  S  V08 

R  S  U 

FPONE 

C  S  V  02 

0  10 

139 

0  317 

C  S  V06 

FSB 

IIOPI 

Oil 

14  0 

00  4  6 

NZE 

C  S  V  0  4 

0  12 

141 

0  17  5 

BM  1 

C  S  V  06 

0  13 

14  2 

017  8 

F  A  0 

0  N  E  P  1 

0  14 

14  3 

0  4  3  1 

NZE 
S  T  U 

C  S  VO  7 

C_S  V  0  H 

01  S 
0  16 

14  4 
14  5 

0  313 

cVvo-T 

C  S  V  0  9 

00  16 

R  SU 

FPONE 

0  17 

14  6 

00  3  3 

S  T  U 

C  S  V  10 

0  18 

147 

0  19  5 

S  T  U 

C  S  V  11 

0  19 

14  8 

045  3 

STL 

C  S  V  1  2 

C  S  V  13 

0  20 

14  9 

00  6  1 

EJ  0  C  A 

S  T  0. 
NZE 

C  S  V.02 

B-EG 

IN   SIN 

0  2  1 

150 
151 

018  3 

C  S  V  02 

0189 

BM  1 

C  S  V14 

C  S  V  1  S 

02  2 

152 

009  2 

C  S  V  1  4 

FAD 

IIOPI 

023 

153 

03  4  5 

NZE 

C  S  V  02 

0  2  4 

154 

032  5 

BM  1 

C  S  V  14 

0  25 

155 

0  228 

.  .    

FSB 

ON  EP  1 

0  26 

156 

0  38  2 

NZE 

C  S  V  16 

C  S  V  02 

0  2  7 

15  7 

026  3 

C  S  VI  5 

FSB 

IIOPI 

02  8 

158 

009  6 

NZE 

C  S  V  02 

029 

159 

027  5 

BM  1 

C  S  V  1  5 

0  30 

160 

0278 

FAD 

0  N  E  P  1 

0  31 

16  1 

0  4  8  1 

NZE 

C  S  V  16 

CSV  02   

.     0  3  2 

162 

0  313 

CSV16 

S  T  U 

C  S  V  0  9 

0  33 

163 

00  6  6 

R  SU 

8  003 

0  34 

164 

007  3 

S  T  U 

C  S  V  10 

0  3  5 

165 

0  5  3  1 

ST  U 

C  S  V  11 

0  36 

166 

0  SO  3 

LO  D 

FPONE 

037 

167 

0  111 

S  TO 

C  S  V  13 

CS  V  13 

0  3  8 

168 

0  19  3 

C  S  V  13 

R  A  U 

C  S  V  12 

0  39 

169 

00  18 

F  A  0 

FPONE 
C  S  V  17 

0  40 
0  41 

170 
171 

00  19 

S  T  U 

026  7 

FAD 

FPONE 

0  4  2 

17  2 
17  3 

0  32  5 

S  T  U 

C  S  V  12 

0  43 

0317 

R  SU 

C  S  V  10 

0  4  4 

174 

0068 

FM  P 

C  S  V  09 

0  4  5 

175 

0  30  5 

F  MP 

C  S  V09 

0  4  6 
0  47 

176 
177 

0070 

FO  V 

C  S  V  17 

0  120 

F  0  V 

C  S  V  1  2 

0  4  6 

17  8 

0  3  2  3 

STU 

C  S  V  10 

0  49 

179 

0  115 

_ 

FO  V 

C  S  V  11 

0  50 

180 

0553 

RAM 

8  003 

0  51 

18  1 

0  15  0 

R  A  U 

8  0  02 

0  5  2 

182 

0  16  5 

FSB 

SIZES 

0  53 

183 

0  12  3 

BM  1 

C  S  V  18 

0  5  4 

184 

0  60  3 

R  A  U 

C  S  V  1  1 

0  5  5 

185 

00  0  7 

F  AO, 
STU 

C  S  V  10 

0  5  6 
0  5  7 

186 
187 

036  3 

C  S  V  1  1 

C  S  II  i  3 

0  328 

C  S  V  18 

R  A  U 

C  S  V  1  1 

C  3  V  02 

05  8 

188 

00  5  6 

F  P  ONE 

10 

0000 

00  5  1 

059 

189 

00  4  0 

F  P  T*  0 

20 

00  00 

00  5  1 

060 
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RESULTS  AKD  DISCUSSION" 
Steady-state  Tubular  Plow  Reactor 

The  concentration  profiles  for  the  various  combinations  of 
M  and  R  for  the  first-order  chemical  reactions  were  calculated  by 
use  of  equations  (9)  through  (12),  (16)  and  (18).   Typical  results 
corresponding  to  R=l/2  for  B.C. (I),  (II)  and  (III)  are  shown  in 
Fig.  1  and  those  corresponding  to  R=2  for  B.C. (IV)  are  shown  in 
Pig.  2.   The  figures  show  that  the  concentration  profiles  obtain- 
ed by  use  of  B.C. (II),  (III)  and  (IV)  were  inconsistent  with  the 
limits  of  both  a  plug  flow  reactor  and  a  uniformly  stirred  tank 
reactor  for  all  values  of  M  and  R;  whereas  the  profiles  computed 
by  use  of  B.C. (I)  were  consistent  with  such  limits. 

The  word  "inconsistent"  means  that  a  computed  conversion  or 
concentration  profile  either  less  than  that  of  a  uniformly  stirred 
tank  reactor  or  greater  than  that  of  a  plug  flow  reactor  is  un- 
reasonable.  Models  which  predict  this  are  therefore  "inconsistent". 
The  numerical  values  of  the  limits  of  inconsistency  were  obviously 
calculated  by  use  of  equations  (16)  and  (18)  for  a  uniformly  stirred 
tank  reactor  and  a  plug  flow  reactor  respectively. 

Figure  3,  however,  indicates  that  the  over-all  conversions 
computed  by  use  of  the  different  sets  of  boundary  conditions  were 
numerically  in  agreement  for  large  values  of  M  (small  degree  of 
axial  dispersion)  and  deviate  increasingly  from  each  other  as  the 
degree  of  axial  dispersion  increases.  Again  the  overall  conver- 
sions computed  by  use  of  B.C. (II)  and  (III)  gave  values  lower  than 
those  obtained  for  a  uniformly  stirred  tank  reactor,  and  those 
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FIG.  I.   Concentration   profile  of  react- 
ant  for  various  M  for  first-order  reac- 
tion with  R*tt  . 
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as  a  function  of  R  for  first-order  reaction. 
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computed  by  use  of  B.C. (IV)  gave  values  higher  than  those  obtained 
in  a  plug  flow  reactor,  especially  for  low  conversions  and  a  high 
degree  of  axial  dispersion,  i.  e.,  low  values  of  R  and  M0   Reactors 
falling  in  this  category  are  sometimes  met  in  laboratory  and  pilot 
plant  experimentation,  though  not  commonly  required  for  design  pur- 
poses. 

The  inconsistency  shown  in  Figs.  1  and  2  led  to  the  establish- 
ment of  the  domain  for  which  any  combinations  of  M  and  R  are  incon- 
sistent, with  the  limits  of  the  domain,  as  shown  in  Figs,  l\.   and  5» 
The  curves  in  Fig.  h   show  the  locus  representing? 

(fii:i)>  =1  "  fM=0  ™  ° 
The  domain  below  the  curves,  then,  corresponds  to  the  condition; 

(fri)  ^^  -  rM^0   <  0 

(fIIl}  ^~1  "  fM=0  <  ° 
The  curve  in  Fig.  5>  indicates  the  locus  representing? 

(fiv)  7=1  "  %=a>  s    ° 
The  domain  below  the  curve,  then,  corresponds  to  the  condition; 

Where  f  represents  the  overall  conversions. 

The  inconsistency  of  the  overall  conversions  using  B.C. (IT) 
has  been  indicated  by  Wehner  and  Wilhelm  (2[|_) .   However,  it  has 
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never  been  established  quantitatively  as  a  function  of  the  system 
characteristics 5  1.  e.,  M  and  R.   The  cases  of  B0Co(III)  and  (IV) 
have  been  hitherto  unknown.   The  overall  conversions  computed  by 
use  of  B.C. (II)  and  (III)  are  consistent  with  the  upper  limit  of 
a  plug  flow  reactor,,  and  those  by  use  of  B.C.  (IV)  are  consistent 
with  the  lower  limit  of  a  stirred  tank  reactor. 

The  design  of  an  isothermal  reactor  should  not  be  based  only 
on  the  overall  conversion  when  the  reactions  to  be  carried  out  in 
the  reactor  are  exothermic  or  endorthermlc.   The  estimation  of  the 
heat  transfer  area  to  be  provided  requires  a  knowledge  of  the  con- 
centration profile,  or  the  conversion  along  the  length  of  reactor. 
Figs,  6  and  7  show  the  comparison  of  concentration  profiles  based 
on  the  different  sets  of  boundary  conditions  at  both  high  and  low 
degree  of  dispersions.   The  figures  supplement  the  previous  state- 
ment that  the  deviations  among  the  profiles  decrease  as  M  increases 
This  means  that  the  effect  of  boundary  conditions  become  negligible 
for  high  values  of  M,  say  of  the  order  of  15,  and  the  reactor  can 
be  assumed  to  behave  as  a  plug  flow  reactor  with  reasonable  ac- 
curacy.  The  percentage  deviations  of  the  overall  conversions  com- 
puted by  use  of  B.C. (II),  (TIT)  and  (IV)  from  those  computed  by 
use  of  BoG.(I)  are  shown  in  Figs,  8,  9  and  10.   The  reason  that 
B.C.  (I)  is  used  as  a  basis  for  the  computation  is  that.,  as  point- 
ed out  previously,  this  yields  conversions  which  are  consistent 
with  both  limits  of  the  reactor.   These  figures  show  that  the  de- 
viations increase  as  M  and  R  decrease,  and  B.C. (II)  always  yields 
larger  deviations  than  B.C. (Ill)  and  (IV).   Since  the  solutions 
to  B.C. (Ill)  and  (IV)  are  in  the  form  of  a  simple  exponential 
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Fig.  6.    Comparison    of   concentration   profile 
of  reacting   systems  (R=2)    based    on 
different    boundary    conditions    at 
high   degree    of    dispersion (M  =  I) . 
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Fig.  7.     Comparison    of    concentration    profile 
of   reacting   systems  (R=2)    based    on 
different    boundary   conditions    at 
low    degree    of    dispersion (M  =  15). 
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Fig.  8.     Percentage     deviations     of     final     conversions 

computed     by     use     of     boundary  conditions   (I) 

from     those     computed     by     use  of     boundary 

conditions  (I)      as     functions     of  M      for     first- 
order     reactions. 
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Fig.    10.     Percentage      deviations     of     final 
conversions     computed     by     use    of     boundary 
conditions  (EZ)     from     those     computed    by     use 
of    boundary     conditions  (I)     as     functions    of 
M     for     first- order    reaction. 
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function9  their  use  to  approximate  B0C0(I)  is  always  more  conveni- 
ent,, and  also  yields  more  accurate  results  than  the  use  of  B0C0(II)o 

The  establishment  of  regions  of  inconsistency  for  the  use  of 
different  sets  of  boundary  conditions  gives  a  set  of  criteria  for 
the  use  of  B0Co  (II) 9    (III)  and  (IV) 0   These  boundary  conditions 
should  be  avoided  in  the  domain  of  inconsistency  for  computing  fi- 
nal conversions o  Due  to  ease  of  solving  simple  exponential  func- 
tions, the  use  of  B0C0(III)  or  (IV)  ,  if  their  limits  are  recognizeds 
offers  the  most  convenient  means  for  computing  the  final  conversions 
The  recognition  of  the  limit s    of  course ,  refers  to  the  domain  above 
the  curves  in  Pigs0  k   and  5<>   The  use  of  B0C0(II)  for  a  homogenous 
tubular  flow  reactor  with  axial  dispersion  is  not  desirable  from 
the  view  point  of  both  accuracy  and  convenience 0   Though  the  use 
of  BoCo(IV)  has  never  been  attempted  previous ly,  it  is  useful  in 
the  computation  of  conversions  around  the  inlet  and  center  of  the 
reactors  since  as  shown  in  Pigse  6  and  7S  the  concentration  profile 
computed  by  use  of  B0C0(IV)  follows  closely  that  computed  by  the 
use  of  BoCo(I)p  except  at  the  exit  of  reactorc   As  stated  previous- 
ly 5  the  simple  analytical  solution  in  exponential  form  for  B0C0(IV) 
provides  a  reasonably  accurate  and  convenient  method  of  evaluating 
con vers ionso 

Unsteady-state  Tubular  Plow  Reactor 

The  numerical  calculations  of  the  impulse  response  function 
tabulated  in  Table  l\.   are  relatively  easy  in  comparison  with  the 
direct  numerical  solution  of  the  original  partial  differential 
equations  equation  (2)0   Some  computational  results  of  the  impulse 
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response  functionsp  both  with  and  without  chemical  reaction,  are 
presented  in  Figs.  11  through  18 «   These  figures  definitely  indi- 
cate the  differences  of  response  based  on  different  sets  of  bound- 
ary conditions.  Pour  sets  of  impulse  response  curves  at  R=0  for 
high  and  low  degree  of  dispersion  are  replotted  in  Pigs,  19  and 
20  to  show  that  as  the  degree  of  dispersion  becomes  small,  say  of 
the  order  of  M=l55  all  the  response  curves  representing  the  four 
boundary  conditions  approximately  converge  to  one  and  the  devia- 
tions among  the  four  response  curves  increase  as  degree  of  dis- 
persion increases o   This  fact  supplements  the  previous  statement 
that  the  effect  of  boundary  conditions  becomes  negligible  above 
the  value  of  M~l5<>   Exactly  the  same  observations  are  given  in 
Pigs o  15  through  18  where  the  response  curves  of  first-order  chem- 
ical reaction  are  considered,  except  that  the  magnitude  of  the  re- 
sponse curves  decrease  rapidly  with  increase  in  R.   This  is  due 
to  the  fact  that  the  magnitude  of  R  is  proportional  to  the  rate 
of  reactions  thus  the  faster  the  reaction  the  faster  the  concen- 
tration decreases. 

One  application  of  such  computed  response  carves  is  the  de- 
termination of  the  axial  dispersion  coefficient  from  a  dynamic  ex- 
periment. A  survey  of  published  literature  shows  that  most  of  the 
methods  used  to  determine  the  mean  dispersion  coefficient  from 
transient  response  experiments  makes  use  of  only  a  single  point  of 
an  entire  experimental  response  curve 9    such  as  the  maximum  point 
or  the  slope  at  unit  time  (Ij_)  (25).   It  is  known  that  the  use  of 
a  single  datum  in  drawing  conclusions  for  an  entire  experiment 
may  lead  to  serious  error  because  of  probable  experimental  error. 
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Fig.  II.      Calculated     response     curves     of      exit 
concentration     against     impulse      for     non-reacting 
systems  (R=0)      based      on      B.  C.  (I)      with     M 

as      parameter. 
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Dimensionless      time      variable,     X 

Fig.  12.    Calculated     response    curves     of    exit 
concentration     against     impulse      for  .  non  reacting 
systems  (R  =  0)      based      on      B.   C.  \H)      with 

M     as     parameter. 
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Fig.  13.       Calculated     response     curves     of     exit 
concentration      against       impulse      for    non-reacting 
systems   (R  =  0)      based      on      B.  C.  (IE)      with 

M     as    parameter. 
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Fig.  14.  Calculated  response  curves  of  exit 
concentration  against  impulse  for  non- reacting 
systems  (R=0)     based    on     B.  C.CE)     with      M 

as     parameter. 
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Fig.  15.     Calculated     response     curves    of     exit 
concentration     against     impulse     for     reacting 
systems  (R=0. 1)      based     on     B.  C.  (I)     with 

M     as     parameter. 
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Fig.  16.      Calculated      response    curves     of 
exit    concentration     against     impulse     for 
reacting     systems  (R=IO)     based     on     B.C. 
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Fig.  17.  Calculated  response  curves  of  exit 
concentration  against  impulse  for  reacting 
systems  (R  =  0.1)     based     on     B.C.  (E)     with    M 

as    parameter. 
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Fig.  19.     Comparisons     of     response     to  impulse 

for    non  reacting     systems  (R  =  0)      based  on 

different     boundary     conditions     at     high  degree 

of     dispersion  (M=l). 
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Another  approach  In  estimating  D  is  reported  by  van  der  Laan  (21), 
Levenspiel  and  Smith  (12)  who  introduced  the  concept  of  the  vari- 
ance of  the  concentration  distribution  by  momentum  analysis ,   This 
method  has  the  advantage  by  avoiding  the  cumbersome  back  transfor- 
mation of  the  Laplace  transform.   The  disadvantage,  however,  is 
that  one  cannot  explain  how  distribution  functions  for  residence 
times  can  be  defind  and  measured  for  actual  systems.   Since  the 
overall  flow  behavior  problem  is  often  of  prime  concern  in  the  de- 
sign of  continuous  processing  equipment,  one  should  compare  response 
curves  calculated  from  the  assumed  models  and  boundary  conditions 
with  the  experimental  curves. 

Since  any  model  of  fluid  dispersions,  including  the  diffu- 
sion-type model,  is  approximate  in  nature,  the  best  procedure  for 
estimating  the  mean  dispersion  coefficient  based  on  any  one  of  the 
models  should  be  the  comparison  of  the  calculated  response  curves 
and  an  experimental  response  curve.   Such  a  procedure  is  illu- 
strated in  Pig,  21, 

The  experimental  data  used  in  Fig,  21  were  obtained  from 
Run  19U  of  Fowler  and  Brown  (6),   Levenspiel  and  Smith  (12)  ob- 
tained the  value  of  mean  dispersion  coefficient  of  2528  ft^/sec, 
for  this  experiment  by  the  methods  of  both  maximum  point  and  vari- 
ance of  concentration  distribution  curve.   Inspection  of  Fig,  21, 
however,  shows  that  the  comparison  of  the  entire  response  curves 
yields  mean  dispersion  coefficient  of  2,7  ft2/sec,  if  the  response 
curve  based  on  B»C.(III)  is  used.   It  also  indicates  that  the  sys- 
tem represented  by  the  experimental  data  may  not  be  represented 
by  B,C,(II)  and  (IV),   In  other  words,  the  comparison  of  the 
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experimental  and  calculated  response  curves  should  lead  to  the  es- 
tablishment of  the  appropriate  models  or  boundary  conditions 0 

There  are  no  data  available  on  impulse  experiments  accompanied 
by  first-order  chemical  reaction.   If  the  data  were  available,  the 
calculated  response  curves  based  on  first-order  chemical  reaction 
(Figs,  15  through  18)  could  be  useful  to  extract  transport  para- 
meters from  a  dynamic  experiment. 

The  typical  results  for  step  input  response  curves  are  shown 
in  Figs.  22  through  25  both  with  and  without  chemical  reaction* 
Due  to  the  consistency  observed  in  the  steady-state  concentration 
profile  with  the  limits  of  both  stirred  tank  and  plug  flow  reactor, 
only  the  response  curves  based  on  B,C„(I)  are  presented  here.   The 
response  curves  with  first-order  chemical  reaction  do  not  approach 
one  as  "£  approaches  infinity  because  of  the  steady-state  conver- 
sion, i.e0<,  the  response  curves  approach  steady-state  conversion 
as  't  approaches  infinity0 

One  method  of  determining  D  from  step  input  response  curves 
is  illustrated  by  Levenspiel  and  Smith  (12)6   They  replotted  the 
data  of  Von  Rosenberg  (23)  on  probability  paper  to  obtain  a  straight 
line.   From  the  property  of  normal  distribution  that  the  interval 
between  the  16th  and  81|th  percentile  points  is  two  standard  devi- 
ations ,  they  evaluated  the  mean  dispersion  coefficient,,   This 
method  requires  an  extra  step,  and  also  involves  appreciable  error 
because  the  actual  experimental  data  does  not  give  a  symmetric  nor- 
mal distribution  curve,   The  best  procedure  for  determining  D, 
again  therefore,  should  be  by  matching  the  calculated  response 
curves  and  an  experimental  response  curve. 
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The  results  of  frequency  response  curves  based  on  B.C. (I) 
(equation  31)  are  shown  In  Figs „  26  through  28  Ln  the  form  of 
dimensionless  Bode  plots.   Frequency  responses  of  the  two  limit- 
ing conditions  of  the  dispersion  model  based  on  equations  (38)  and 
(1|3)  are  also  shown  in  the  figures  for  the  sake  of  comparison. 
Only  the  frequency  response  curves  based  on  BoC.(I)  are  presented 
herewith  because  of  the  consistency  observed  in  the  previous  sec- 
tion with  the  two  limiting  cases.   Other  frequency  curves  based 
on  different  sets  of  boundary  conditions  are  presented  only  for 
the  sake  of  comparison  in  Figs.  29  through  3^« 

Figure  26  shows  that  the  plots  of  amplitude  ratios  and  phase 
shifts  change  monotonically  between  the  two  limits  for  the  nonre- 
acting  systems  in  the  range  of  6$Q     which  is  of  practical  Interest 
(15)«   This,  however,  is  no  longer  true  for  the  reacting  systems, 
as  shown  in  Figs 0  27  and  28.   As  the  relative  reaction  rates  char- 
acterized by  R  increase,  amplitude  ratio  decreases  rapidly,  es- 
pecially at  large  values  of  M.   For  example,  at  the  conditions  for 
which  R  and  M  are  greater  than  6  and  10  respectively,  the  ampli- 
tude ratio  is  very  small  and  practically  is  negligible 4   Under 
these  conditions  of  rapid  reaction  and  small  degree  of  dispersion, 
a  cyclic  unsteady-state  system  has  little  effect  on  the  perform- 
ance of  the  system;  and  therefore,  the  frequency  response  data  are 
of  little  use.   In  contrast,  the  phase  shift  data  are  only  slight- 
ly affected  by  the  change  in  the  values  of  M  and  R.   Thus,  it  can 
be  said  that  the  phase  shift  data  are  of  little  value  in  deter- 
mining the  characteristics  of  the  system. 

The  performance  function  represented  by  equation  (118)  (See 
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Fig.  26.     Frequency     response     diagram     of 

non-reacting     system    (R  =  0)      based 

on      B.C.  (I)     with      M      as      parameter. 
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Fig.  27.     Frequency      response      diagram      of 
reacting      system     (R=2)      based      on 
B.    C.    (I)      with      M      as      parameter. 
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Appendix  IV)  was  originally  derived  by  Kramers  and  Alberda  (10) 
for  non-reacting  systems 0   In  order  to  avoid  complexities  in  trans- 
forming equation  (118)  to  the  computable  form  of  equation  (31) ,  it 
was  suggested  that  equation  (31)  be  approximated  by  the  perform- 
ance function  of  n  numbers  of  the  perfectly  mixed  vessels  connected 
in  series o   Such  approximation  is  good  only  for  large  values  of 
M(c0a0  M=  15) o   Very  frequently  approximate  boundary  conditions 
which  do  not  conserve  the  continuity  requirements ,  i0e0,  BoCe(II) 
and  (III)  j  are  employed,,   Kramers  and  his  coworkers  (11)  employed 
BoCo(II)  and  (III),  and  the  evaluation  of  the  performance  function 
based  on  B0C0(III)  (equation  108  of  Appendix  IV)  has  been  frequent- 
ly presented  in  the  literature  (7)  (22) 0 

The  amplitude  ratios  and  phase  shift  based  on  these  different 
sets  of  boundary  conditions  are  compared  in  Figs„  29  through  3h 
both  with  and  without  chemical  reaction  at  high  and  low  degree  of 
dispersion,.   As  can  be  seen  from  Pig.  29 9    the  agreement  of  the 
amplitude  ratios  and  phase  shift  curves  for  non-reacting  systems 
based  on  approximate  boundary  conditions  are  quite  poor  at  low 
values  of  M  (large  degree  of  dispersion) 0  An  increasing  degree 
of  agreement,  however,  is  observed  in  Pigo  30  as  M  increases „   The 
deviations  among  the  phase  shift  curves  are  very  small  when  M  is 
approximately  10  or  greater,  and  the  amplitude  ratio  curves  prac- 
tically coincide  when  M  is  approximately  15  or  greater,,   It  is 
also  seen  that  phase  shift  curves  are  generally  less  affected  than 
amplitude  ratio  curves.   The  performance  function  curves  of  the 
reacting  system,,  shown  In  Figs0  31  through  3U-.,  reveal  the  effects 
of  using  different  sets  of  boundary  conditions  more  clearly  than 
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Fig.  29.     Comparison     of     frequency     response 
diagrams     of     non-reacting     system  (R  =  0)      based    on 
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boundary     conditions     at     high     degree     of     dispersion, 
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the  corresponding  curves  for  the  non-reacting  systems  ~-  the 
greater  the  value  of  R,  the  greater  the  effect  of  the  boundary 
conditions.   It  is  also  seen  from  these  figures  that,  at  any 
given  relative  reaction  rate  characterized  by  R,  the  effect  of 
different  sets  of  boundary  conditions  becomes  smaller  as  the  de- 
gree of  mixing  decreases 0 

An  illustration  of  the  determination  of  the  mean  dispersion 
coefficient  from  a  frequency  response  experiment  for  non-reacting 
systems  by  comparing  available  experimental  data  (7)  (10)  (22) 
with  calculated  values  is  presented  in  Pigs.  35>  through  37»  These 
figures  show  that  the  use  of  approximate  boundary  conditions  can 
give  rise  to  an  appreciable  error  in  determining  the  degree  of 
dispersion,  especially  at  low  values  of  M.   As  has  been  discussed 
in  the  previous  section,  the  use  of  B.C. (Ill)  is  quite  desirable 
under  some  circumstances  due  to  the  simplicity  of  equation  (33) o 
For  the  preliminary  design  of  a  tubular  flow  system  and  the  syn- 
thesis of  a  complex  control  system  in  which  the  use  of  a  compli- 
cated performance  function  is  involved,  the  use  of  approximate 
boundary  conditions  should  be  avoided. 
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CONCLUSION 

The  general  extent  of  axial  dispersion  and  its  effect  on  the 
performance  of  isothermal  tubular  flow  reactors  can  be  character- 
ized  by  a  diffusion -type  model  with  appropriate  boundary  conditions, 
Among  the  four  boundary  conditions  considered,  only  B0C0(I)  yields 
conversions  and  concentration  profiles  consistent  with  the  theo- 
retical limits  of  both  a  uniformly  stirred  tank  reactor  and  a  plug 
flow  type  reactor,  and  for  this  reason,  the  use  of  BoC0(II),  (III) 
and  (IV)  should  be  avoided  in  and  around  the  domain  of  inconsistent 
conversions  based  on  a  one -dimensional  transport  equation,, 

There  are  no  data  available  on  transient  response  experiments 
accompanied  by  first-order  chemical  reaction*   If  the  data  were 
available,  the  calculated  response  curves  based  on  first-order 
chemical  reaction  could  be  used  to  determine  the  dispersion  co- 
efficient by  the  same  method  presented  in  this  research  for  the 
non-reacting  case0 

The  response  curves  presented  in  this  research  for  various 
types  of  transient  input  should  be  useful  for  practicing  process 
control  engineers  and  for  those  who  wish  to  extract  the  transport 
parameters  from  dynamic  response  experiments. 
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NOMENCLATURE 


95 


(A.R.) 


=  angle  defined  as  J-  Tan"-1- 


2(  CdS  ) 

2R  +  M 


,  degree. 


=  a  constant  in  equation  (27)  which  is  numerically 
equal  to  the  steady-state  conversion,  dimension- 
less. 

=  amplitude  ratio,  dimensionless. 


1  +  2R)2  +  k(  <^)2 

M         M2 


B 
C 

cx(t) 
c2(t) 

D 
D 

*l 

fM=0 

3 

k 
L 
M 
n 
P.P. 


a  constant  in  equation  (28)  which  is  numerically 
equal  to  the  steady-state  conversion,  dimension- 
less. 

concentration  of  reactant,  m/1. 

inlet  concentration  of  reactant,  m/1. 

exit  concentration  of  reactant,  m/1. 

axial  dispersion  coefficient,  l2/t. 

mean  value  of  D,  l2/t. 

conversion  calculated  based  on  B.C.(i) ,  dimen- 
sionless. 

conversion  calculated  based  on  a  stirred  tank 
reactor,  dimensionless . 

conversion  calculated  based  on  a  plug  flow 
reactor,  dimensionless. 

reaction  rate  constant,  t""^. 

length  of  reactor,  1. 

uL/2D,  1/2  axial  of  Peclet  number,  dimensionless 

order  of  reaction. 

performance  function  defined  as  the  ratio  of 
Fourier  transform  of  output  to  input. 
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Q  =  amount  of  tracer  Injected  In  impulse  function,  m. 

R  =  kL/u,  reaction  number,  dLmensionless. 

r(C)  =  rate  of  reaction. 

s  =  function  of  Laplace  transform, 

t  =  time. 

u  =  average  flow  velocity  of  reacting  mixture,  l/t. 

v  =  volumetric  flow  rate,  l3/t„ 

V  =  volume  of  reactor,  l30 

W(s)  =  transfer  function  defined  as  the  Laplace  trans- 
form of  output  to  input. 

x  =  Mb  cos  a 

y  =  Mb  sin  a 

z  =  co-ordinate    in  the   direction   of  flow,    1. 

Greek  Letters 


o^-,  =     sinh  x   cos  y   -    ^  ^P  '    cosh  x  siny 


7T 
+  b(cos   a   cosh  x  cos  y   -   sin  a   sinh  x  sin  y) . 

^2  =     f  )  sinh  x  coy  y  +   cosh  x  sin  y 

+  b(cos   a   sinh  x  sin  y  +  sin  a   cosh  x  cos  y)  e 

c^o  -     sinh  x  cos  y 

+  b(cos   a  cosh  x  cos  y   -   sin  a  sinhx  sin  y) . 

o(|,  =     cosh  x  sin  y 

+  b(sin  a   cosh  x  cos  y  +  cos   a   sinh  x  sin  y) . 


/S  =  (1  +   2R/M)2    . 

rs'  r-  i 

p  =  _1  +  (26/M)(s   +  k)]    2 

(S  =  Ql  +    (29/M)(J60+  k)J  * 
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^  n     =   roots  of  cot  o     -  h   (  —  -  — r=-  ) 

K   M     S  J 

<3  „     =  roots  of  tan  d     =  -  — — 
n  M 


=  phase  shift,  degree 


4> 

\  -  z/L,  dimensionless  distance  from  the  inlet  of 
reactor, 

y 

0  -  C/C]_,  dimensionless  concentration. 

^  =  the  dimensionless  time  variable,  \,/q    . 

'Cq  =  thickness  of  impulse  in  time  domain,  dimensionless 

6  =  L/u,  t. 

^0  =  angular  frequency,  rad./t,, 


98 


APPENDIX 
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APPENDIX  I 
Derivations    of  Steady-state   Solutions 

The  solutions  of  equation  (8)  which  satisfy  the  boundary  con- 
ditions (I),  (II),  (III)  and  (IV)  are  obtained  by  means  of  Laplace 
t  rans  f  orma  t  i  ons  <, 

Laplace   transform  of   equation    (8)   with   respect   to  >£    is: 

/(s)    _  s    ^(0+)    +      /(0+)    -   2M    /(0+) 

s2   -  2Ms    -   2MR  (U8) 

For  the  solution  using  B«C,(I),  the  Laplace  transform  of 
equation  (3)  is  substitued  in  equation  (i|8)  to  obtain: 

/(S)  =  s  /(0+)  -  2M 

s2  -  2Ms  -  2MR  (J4.9) 

By  means  of  partial  fractionation,  equation  (l±9)    is  sepa- 
rated into  two  parts: 


yfSN  .  il±jLj   ^(q+)  - 2    d-/g)  /(q+)  -2 

'  2/5[s    -  M(H*y?   )]  2jp    -  M(l-^    q- 

The    inverse    transform   of   equation    (50)    gives: 

Av  =  12*01  f^)  -2  exp[-M(1+/6)n 

d-/5)  y(o+)  -  2 


(5o) 


2 


exp[>(l- /$)>£]  (£l) 


According  to  boundary  condition  represented  by  equation  (1+) , 
differentiation  of  equation  (5>1)  and  setting  it  equal  to  zero 
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gives  : 

/(0+)    =  2{(l+l)exp[M(l-f/3)3   -    (l-/3)expfa(l-/6i]} 
( 1+  j3)2exp[vL{  1+  /S )]  -    (1-  /S)2exp['M(  1-/3  )] 

Substitution  of  equation  (52)  into  (51)  and  simplification 
gives  equation  (9)  of  text. 

For  the  solution  using  B.C. (II),  the  Laplace  transform  of 
equation  (5)  is  substituted  into  equation  (I48)  and  simplified  to 


/(8)  =  (s  -  2M)  +  /(0+) 


75 (53) 

3d   -  2Ms  -  2MR 


By  means  of  partial  fractionation,  equation  (53)  is  sepa- 
rated into  two  parts  as: 


J.    >  ._  -  M(l-l)  +  y'(0+)         -M(l+/S)  +  /(0+) 
'       2M/S  [s-M(l+/S  )]       2M^[s-M(l-/5)] 


(A) 


The  inverse  transform  of  equation  (5^)  yields: 

According  to  boundary  condition  represented  by  equation  (I;), 
equation  {^)    is  differentiated  in  order  to  solve  for  /(0+)„ 
It  gives: 


/(0+)    =  M(l-^2) 


exp[M(  1+  £  )]   -   exp[k(  1-/6)3 


(l+/3)expCM(l+/3  fl-    (l-/6)exPrM(l-/3)]J 
r  ^  r      (56) 
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Substitution  of  equation  (56)  Into  (55)  and  simplifica- 
tion gives  equation  (10)  of  text. 

For  the  solution  using  BoC0(III),  substitution  of  equa- 
tion (5)  into  equation  (lj.8)  gives  equation  (55)  <,     As  )7 
approaches   infinity  in  equation  (^5) ,    the  second  term  in  the 
right-hand  side  of  equation  {S^>)    vanishes  as? 

M(  &   -  1)  ±      /(0+)  (5?) 

2M£  °M 

Solving  for  J  (0+)  In  equation  (57)  and  substituting  it  in- 
to equation  (55)  *  equation  (11)  of  text  is  obtained. 

For  the  solution  using  B.C. (IV),  substitution  of  equation 
(3)  into  equation  (1^8)  gives  equation  (51).   As  Yl     approaches 
infinity  in  equation  (51),  the  first  term  in  the  right-hand  side 
of  equation  {%!)    vanishes,  hence 

0  -  (^  /t°+>  -2  (58) 

2/3 

Solving  for  /(0+)    in  equation  (58)  and  substituting  it 
into  equation  (51)  yields  equation  (12)  of  text. 

APPENDIX  II 

Derivations  of  the  Transfer  Functions,  W(s) 

The  Laplace  transform  of  equation  (2)  with  respect  to  t  at 
zero  initial  concentration  is: 

aC(z,8)  =  D-2f2i|ial  -  tU£l£^i-  kc(z,s)  (59) 

£  z^  d  z 
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Equation  (59)  is  in  the  form  of  a  second-order  differential 
equation  whose  general  solution  is: 


§  (Ufi)  f  (1-/3) 

C(z,s)  =  Aie  dD  +   A2e  ^  (60) 


where  An,  A~  =  arbitrary  constants. 

The  transfer  functions,  defined  as  the  ratio  of  the  Laplace 
transforms  of  output  to  input,  then  become 

W(s)  =  C(L?s)  (61) 

C1(0,s) 

In  order  to  solve  for  the  arbitrary  constants  in  equation 
(60),  the  boundary  conditions  tabulated  in  Table  2  are  trans- 
formed as: 


B.C. (I) 


uC(0,s)  =  uC(0+,s)  -  D  dC(0+,s)  (62a) 

dz 

aC'L|3)  =  0  (62b) 

o  z 


B.C. (II) 


uC(0,s)  a  uC(0+s)  (63a) 

(L,i 


ac(L's)  =  0  (63b) 


uC(0,s)  =  uC(0+,s)  (61j.a) 

B.C. (Ill) 

lim  C(z,s)  =  0  (6U-'b) 


B.C. (IV)  ' 


fuC(0,s)  =  uC(0+,s)  -  D  3C(0+,s)  (65a) 

B  z 

lim  C(t,s)  =  0  (65b) 

Z  -»-co 
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For  the  transfer  function  corresponding  to  B.C. (I),  equation 
(60)  is  differentiated  with  respect  to  z: 


uz 


3C(z,s)  _  Aiu   l<1+0\  A£  e  §  (1-/3') 

2D 


az       2D 
where   £'  =  Q.  +  (^-)(s  +  k)] 


(66) 


According  to  equation  (62b),  equation  (66)  is  set  to  zero 
to  solve  for  one  constant  in  terms  of  another  as: 


An  =  - 


(67) 


where  M  =  H£ 
2D 


Using  equation  (62a),  the  Laplace  transform  of  the  input  is 
obtained  as: 


C!(0,s)  =£1(1-^)  +  t£il+p') 


(68) 


Substitution  of  equation  (67)  into  equation  (68)  gives: 


C^CC^s)  = 


A2 
2~ 


,    M(l+/S)       ,  _     Md-^T 
(1+0)2  e       ,  (1,^)2  e    r 


(1+/T)  eM<1  +  /^ 


(69) 


The   Laplace    transform  of   output   is: 


C(L,s)    a  A]_e 


M(l+ 16  ) 


+  A2© 


M(l-/5) 


(70) 
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Substitution  of  equation  (67)  Into  equation  (70)  gives: 


C(L,s)  =  A2- 


y   M(l-/) 
+  2P  e 


1  +/3/ 


fi' 


(71) 


According  to  equation  (61),  the  transfer  function  is  the 
ratio  of  equations  (71)  to  (67),  which  upon  simplification  gives 
equation  (19)  of  text. 

For  the  transfer  function  corresponding  to  B.C. (II),  equa- 
tion (63a)  is  substituted  into  equation  (60)  to  obtain: 


C;l(0,s)  =  A]_  +  A2 


(72) 


Since  the  second  condition  of  B.C. (II)  is  the  same  as  that 
of  B.C. (I),  equation  (67)  still  holds  for  B.C. (II),  and  substi- 
tution of  equation  (67)  into  (72)  gives: 


Cx(0,s)  =  A2 


,   M(l+/S)        ,       M(l-/6) 
(l+fi)    e        -  (1-/5)  e 

~ — m(  i+  #j 


(1+/6  ) 


(73) 


From  equation  (61),  the  transfer  function  corresponding  to 
B.C. (II)  is  the  ratio  of  equations  (71)  to  (73),  which  upon  sim- 
plification gives  equation  (20)  of  text. 

For  the  transfer  function  corresponding  to  B.C. (Ill), 
equation  (61j.b)  is  substituted  to  equation  (60)  to  obtain: 


lim  C(z,s)  =  0  =  Ax(co)  +  A2(0) 


(7k*) 


"OO 


io5 


So  An  mast  be  zero  In  order  to  satisfy  equation  (6I|-b). 
Ai  =  0  (7i;b) 

Since  the  first  set  of  B.C. (Ill)  is  the  same  as  that  of 
B.C. (II),  equation  (72)  still  holds  for  B.C. (Ill) .   Substitution 
of  equation  (7k)    into  (72)  becomes: 

Cx(0,s)  =  A2  (75) 

The  Laplace  transform  of  the  output  is: 

M(l-/S') 
C(L,s)  =  A2  e  (76) 

The  transfer  function  corresponding  to  B.C. (Ill)  is  the 
ratio  of  equations  (76)  to  (75>)»  which  upon  simplification, 
gives  equation  (21)  of  the  text. 

Since  B.C. (IV)  is  the  combination  of  the  first  set  of  B.C. 
(I)  and  the  second  set  of  B.C. (II),  equations  (68)  and  (7^)  still 
hold  for  this  case.   Substitution  of  equation  (7^-)  into  (68) 
gives  : 

Cl(0,s)=i£ii^2  C77) 

2 

The  transfer  function  corresponding  to  B.C. (IV)  is  the 
ratio  of  equations  (76)  to  (77) »  which  upon  simplification, 
gives  equation  (22)  of  text. 
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APPENDIX  III 

Systematic  Derivations  of  The  Response  Functions  of  The  Exit  Con- 
centrations Against  Impulse  Function  and  Step  Inputs. 

For  the  response  functions  of  the  exit  concentration  against 
impulse,  it  is  known  by  L* Hospital  rule  that 

C(L,s)  =  C1  W(s)  (78) 

For  the  solution  corresponding  to  B.C. (I),  equation  (19)  is 
substituted  into  equation  (78)  to  obtain: 

C(L,s)  =  Cx  JkE (79) 

-MU-/37)        ,    -M(l+/S') 


(1  +/)2  e     r   -  (l-^)2  e 


The  inverse  transformation  of  equation  (79)  is  carried  out 
by  means  of  the  Heavlside  partial  fraction  theorem.   Since  equa- 
tion (79)  is  in  the  form  of: 

q  •  ( 5) 

C(L,s)  =  21 (80) 

PL(s) 

where   q.  (s)  =  I4.C  p 

-m(i-^)  ,      -mi+p') 

p  (s)  =  (1+^)2  e  -  (l-/?)2  e     r 

the  inverse  transform  of  equation  (80)  gives: 


CO 


qt(sk)    M 


C(L,t)  =  >  ^  k'    e  K  (81) 


where  sk  is  the  root  of  P]_(s)  =  0 

Solving  for  the  roots  of  Pi(s)  «  0,  one  obtains 
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sk  = 


Sn     +  M 
2M6 


-  k 


where 


<3n  =  root  of  cot  d>     » 


*(■ 


M 


M  \ 

s  J 


(82) 


Differentiation   of  Pt(s)    with  respect   to  s   gives: 

f(s)    =   -  -^-Tm  cos   M*<      +  2   cos   M*<       -   2M^     sin  M<A 

M*    L 


(83) 


-   2<*      sin  Mo<      -  McK  2   cos   M^l 
where         ©(=    jfi 

Substitution  of  equation  (82)  into  the  expressions  of 
P]_(s)  a^d  q]_(s)  gives  respectively: 

M 


q(sv)    =   C- 


P(sk)    = 


<Tr 


2   c>  ne 


M 


(8U) 


En 


Mcos    d  n  +  2  cos     on  -  2    on  sin    on 


2    S"n  _._    r_  ^"n        _      ^ 


M 


sin 


<S"n  - 


COS 


M 


■] 


(85) 


Substitution  of  equations  (8i|.)  and  (85)  into  equation  (81) 
gives : 
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C(L,t)    _^> 1 

Cx  "'mcos  ^"n+2cos  £n-2  £nsln  Sn~  —  sin  Sn- 


cos 


M 


M 


x 


2  <TneM\  (_Jn 

.  m    )  \    e 


exp 


Sn    ,-r  M2 


t   -  kt 


2Me 


<5"n 
(86) 


Multiplication   of  both  numerator  and    denominator   on  the 
right-hand  side    of   equation    (86)    by    (M  sin  Sn  +    Sncos  on)    give; 


C(L»t) 
Cl 


CO 


n=i 


2  on    \     /m  sin  on  +     Sn  c os  Sn 


M 


x  exp 


Si 


+   M£ 


2M9 


K 


t    -  kt 


(87) 


where  K  =    (Mcos  d~n+2cos  <5n-2  cTnsin  Sn-Z-^sin  dn° 

M 

x    (Ms in  <fn+  Jncos  <3n) 


^fl-cos  d  n) 


M 


K 


Simplification   of  K  gives 

=  -  JjL      ( M2  +  S%.     +  2M) 

M 


^/ 


0 


(88) 


Substitution  of  equation  (88)  into  ^(8?)  and  s implication,, 
gives  t 


C2(L,t)0  ^— ■  2  Jn(M  sin  <fn  +    cTn  cos  cfn) 


n=i 


(M2  +    ^Tn2  +   2M) 


x   exp 


(89) 


M  - 


Sn     +  M2 


2M 


+R 


i)€ 
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4. 

where  £  =  — ~ 

The  left-hand  side  of  equation  (89)  is  transformed  in  terms 
of  the  amount  of  tracer  injected,  Q,  and  the  volume  of  reactor, 
V,  by  the  following  relationship: 


rK 


Q  -  u 


C(\)  d\  =  CX1T  d2u  (90) 


v  =  Trail,  (9D 


lrJ 


where  dp  =  diameter  of  reactor „ 


Hence,  substitution  of  equations  (90)  and  (91)  into  the 
left-hand  side  of  equation  (89)  gives  the  following  relation- 
ships 

C2(L,t)      C?(L,t)      L     Cp^d^L   C2V  ,   % 

— =  — —  x  -  «  -^_£-=-£-  (92) 

Cx         Cx        u     C-^dpi     Q 

Using  equation  (92)  in  place  of  equation  (89)  gives  equa- 
tion (23)  of  text. 

For  the  solution  corresponding  to  B.C. (II),  equation  (20) 
is  substituted  into  equation  (78)  to  obtain* 

C(L,s)  =  C1 jJL. (93) 


(l+/8')  e-M(1+£''  -  (l-p)    e-M(1^'> 


Equation  (93)  is  in  the  same  form  as  equation  (80),  except 
p^(s)  and  q-.(s)  are  replaced  by  the  following  two  expressions; 


q2(s)  =  2Ci/S 

p2(s)  =  (!+£')  e     r       -  (!-£)  e 
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Solving  for  the  root  of  p2(s) 
£n  +  M2 


0,  one  obtains: 


Si.  =  - 


-  k 


2M0 


s 


where  <^n  =  root  of  tan  o     =  -  — — 

M 


Ok) 


Differentiation  of  p~  with  respect   to  s  gives: 
p   (s)    =     (-  J|y-)     (M  cos   M<*    +  cos   M*<     -  M*<    sinM<=<) 


(95) 


Substitution   of  equation   (91;)    into  the   expressions    of  P2(s) 
and  q2(s)   gives    respectively: 


*2<*k>=#*M 


(96) 


P2(sk) 


-2-\  (M  cos  a   n  +  cos  <5n  -  <Sn  sin  <5"n)    (97) 
En) 


Substitution  of  equations  (96)  and  (97)  into  equation 
(81)  gives: 


C2(L,t) 


J22_ 


n  =  i 


SI 

M 


exp 


K, 


M  - 


Sn     +  M2 


2M0 


t  -kt 


(98) 


where    K2  =  Mcos<2n+cos<5n-    on  sin  an. 


Simplification   of  K     gives: 


Ill 


K2  = 


M2  + 


M  +  on 


M 


cos 


Sn 


(99) 


Substitution  of  equation  (99)  into  (98),  upon  simplifi- 
cation, gives; 


C2(L,t)e 


CO 


r2 

On 


ftTj     (M2  +  M  +  Sn^)    cos  <Tn 


x  exp 


M 


s 


n  +  M' 
2M 


+  rJC 


(100) 


Substitution  of  equations  (90)  and  (91)  into  the  left-hand 
side  of  equation  (100)  gives  equation  (2lj.)  of  text* 

Pop  the  solution  corresponding  to  B,C0(III),  equation  (21) 
is  substituted  into  equation  (78)  to  obtain? 


C(L,s)  m   Cx  exp[M(l»/6/) 


(101) 


u2 


If  one  lets  g-,  =  L/^B*  and  g2  =  ~=r  +  k  ,  equation  (101)  is 
simplified  to? 


0(L.«)  -  0j  eM  e"Sl  ^ 


+  s 


(102) 


Equation  (102)  is  transformed  to  the  following  expression 
by  the  rule  of  "shift  of  original  function"  as  % 


C(L,t)  m     *'*&  J>   "X  L"*S1  ^ 


(103) 


The  inverse  transformations  of  equation  (103)  is  carried 
out  by  means  of  tables  of  integral  transforms  (5) f 
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A 


-1  (e"Sl  ** 


(101;) 
Substitution  of  gi  and  g2  into  equation  (101;)  gives: 


Op©    /  M       3     r 


(105') 


Use  of  equations  (90)  and  (91)  in  place  of  the  left-hand 
side  of  equation  (10^)  gives  equation  (25>)  of  text0 

For  the  solution  corresponding  to  B.C. (IV),  equation  (22)  Is 
substituted  into  equation  (?8)  to  obtains 

C(L,s)  -  C,   2  exp[M(l-/S  j]  (106) 

U  +  /T) 

,  gj,  «  E-_  +  k,  and  gc^  ^ 

§3 


If  one  lets  g3  a  Ji-  ,  g^  «  2L-.  +  k,  and  g£  -  . 


equation  (106)  becomes 


£ii^i  -  2gs  eM    j  e"83  (^ZZ  t  ( 107) 


Ci  Jg5+ffy 


The  inverse  transformations  of  equation  (107)  by  means  of 
"shifts  of  original  function"  and  the  tables  of  integral  trans- 
forms (5>)  gives t 
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^  -  2%  ^  <k  exp(-  |r)-  s5  -•p(g3«5+s/t) 


x  Erfc 


S3 


2  \Tt 


+  gs-  /t 


(108) 


Substitution  of  g-,  g?  and  g^  Into  equation  (108),  upon  slm« 
pllflcation,  gives  equation  (26)  of  text. 

The  response  functions  of  the  exit  concentration  against 
step  Input  are  known  to  be; 


C(L,s) 


CnW(a) 


(109) 


Since  the  transfer  functions,  W(a),  corresponding  to  each 
of  four  boundary  conditions  and  their  inverse  transforms  are 
known  by  the  derivations  of  response  functions  against  impulse, 
the  following  relationships  can  be  expressed  from  equation  (109) 
by  the  method  of  convolution; 


f 


€ 


C(L,t) 

Cl 


CO 


JQn., 


2  <gn(M  sin  on  4-     Sn  cos  Sp.) 
(M2  +   2M  +    <fn2) 


M- 


(On 


+  lic 


2M 


a 


+  r)X  d\ 


(110) 


(M2  +  M  +  J~n2)  c  os  <Tn 


x  exp 


M 


(JL±j£  +R)x 


(111) 


2M 


d\ 
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C(L9t) 
C, 


C(L,t) 


27T 


xexp   M  -  i  (mA  +  1!  +2Xr)  J  dX 


:|exp[l 


M-(M  +  2R)  Zl|  Erfc 


-  exp(2M  -  R\)  Erf  c  ^ 


M 
2\ 


2A 


d\ 


(112) 


(113) 


The  integration  of  equations  (110)  and  (111)  are  straight 

forward,  and  their  results  can  be  obtained  easily  as  shown  by 

equations  (27)  and  (28)  of  the  text  respectively,,   The  Integra - 

i 
tions  of  equations  (112)  and  ( 113)  are  difficult  to  carry  out 

by  the  method  of  convolution  because  they  contain  the  error  func- 
tions,, Heaviside»s  partial  fraction  theorem  was  found  to  be  more 
convenient  for  the  inverse  transform  of  these  two  equations . 

For  the  solution  corresponding  to  B.C. (Ill),  equation  (21) 
is  substituted  into  equation  (109)  to  obtains 

CieM(lV) 


C(LPs) 


(ilk) 


If  one  lets  g,  =  — ±L-  and  g ._  =  —  +  k,  equation  (lll|.)  is 
simplified  tot 


"86  /£ 


7 


4-  S 


M 


(ll£) 
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Equation  (115)  is  transformed  to  the  following  expression 
by  the  rule  of  "shift  of  original  function"  as: 


M-g?t  »   -1 


e~S6  & 


(116) 


The  inverse  transformation  of  equation  (116)  is  carried  out 
by  means  of  the  table  of  integral  transforms  (5>)  2 


C(L,t)  _  i 
Cl     * 


M-(-g6  /g^) 


Erfc 


( 


B6 

2  (t 


^7*0 


+  e 


w-(-g6  ^) 


(11?) 


Erfc 


(F7?+  ^) 


Substitution  of  g,  and  g   into  equation  (117)  gives  equation 
(29V  of  text. 

The  analytical  solution  corresponding  to  B,C,(IV)  is  found 
to  be  difficult  because  it  contains  the  error  function  and  can 
not  be  integrated  by  the  method  of  convolution  as  shown  in  equa- 
tion (113)o   The  solution  for  a  step  input  corresponding  to 
B.C(IV)  is  left  in  integral  form  as  shown  in  equation  (30)  of 
the  text*   An  approximate  method,  however,  such  as  integration  in 
a  series  form  could  be  carried  out  by  expanding  the  error  func- 
tions „ 
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APPENDIX  IV 

Derivations  of  Frequency  Response  Functions  In  Terms  of  Amplitude 
Ratios  and  Phase  Shift. 


Since  the  performance  functions  are  defined  as  the  ratio  of 
Fourier  transforms  of  output  to  input,  they  are  systematical}/ 
obtained  from  the  transfer  functions  by  replacing  s  by  jW  as? 


P.P.  =  - 


kl 


(1  +  /3")2   expQ-MU-^Tl  -U-/f)2  exp  [-M(l+/d-')] 


(118) 


P.P.  = 


2L 


(1+/S")  expE-M(l-/8VJ  -(1-/5")  exp  [~M(  1+^")] 


P.Fo  -  exp^Md-^'f] 


(119) 


(120) 


P.P.  -  2  e*P  CmU-I*)] 
U+/3") 


where  ^  Jl  +(S^)(jW  +  k) 


(121) 


The  complex  number,^  ,  of  equations  (118)  through  (121)  is 
arranged  in  polar  form  as  2 


L   M       M 


(l+  Sf  4-  kL^Jglf 


Ja 


(122) 


where  a  = 


=  4  Tan-i  (  _2«S_  \ 

>  2R  +  M  / 
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For  the  solution  corresponding  to  B0C0(I),  equation  (122) 
is  substituted  to  equation  (118)  to  obtain-? 

C2(^^>)    m  kP*e      (coaa+Jalna)  (123) 

Cx(j^>)        k\}+  ^~j]sinh(x+jy)    +  i*.b(  cosa+jsina)cosh(x+jy) 

Separation  of  the  denominator  in  equation  (123)  into  a  real  part 
and  an  imaginary  part  gives* 

M 
Co(j*60)   lj.be   (cos  a  +  j  sin  a) 

2      -  . ■ - ( 121+) 

C1(jo3)      l+(  ^  +  j  c<2) 

where  <^_  and  ^p  are  given  in  the  text. 

Equation  (12l+)  is,  then  transformed  into  the  polar  form  as* 


M 


c2(j*»      be 


2\2 


x(j^)    (^  +  o<2  ) 


a-Tan 


(5fL 


(125 


Transformation  of  equation  (125)  into  the  amplitude  ratio 
and  phase  shift  results  in  equation  (31)  of  the  textu 

For  the  solution  using  B0C.(II),  equation  (122)  is  substi- 
tuted into  equation  (119)  which,  upon  simplification,  gives: 


C  (  1^0)  M 

2      _  2be   (cosa  4-  jsina)  H?AX 

CX(JA))  "{l+b(e^a)}  exp(x+jy)  -{l-b(eJa}}  exp[-(x+jy)] 

Further  separation  of  the  denominator  in  equation  (126) 
into  real  and  imaginary  part  gives; 
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M 


C2(j£0)    2beM  (cos  a  +  j  sin  a) 


(127) 


C1(j«0) 


2(  cK3  +  ^) 


where  &(->    and   ^j  are  given  In  the  text. 

Rearrangement  of  equation  (127)  into  the  polar  form  gives: 


C2(j60)       be 


M 


exp 


2x2 


C-lUo))    (  ^   +  *<^) 


a  -  tan" 


1  (^) 


(128) 


Equation  (128)  is  then  transformed  into  amplitude  ratio  and 
phase  shift,  which  give  equation  (32)  of  the  text. 

Substitution  of  equation  (122)  into  equation  (120)  gives  the 
solution  for  B.C. (Ill)  as: 


Cp(j^>) 

*l[l-b( 


»  exp  M  l-b(cos  a  +  j  sin  a 


Cx(j60) 


)]  = 


M~x  -jy 
e   e 


(129) 


Equation  (33)  of  the  text  is  obtained  by  transforming  equa- 
tion (129)  in  terms  of  amplitude  ratio  and  the  phase  shift, 

For  the  solution  using  B.C. (IV),  equation  (122)  is  substitu- 
ted into  equation  (121),  which  upon  simplification,  gives: 


c2(ja))      M-x 

— —  a  2  e 

ox(3a)) 


cos  y  -  j  sin  y 

(l  +  *)+  j(sin  a) 

W 


(130) 


Rearrangement  of  equation  (130)  in  polar  form  gives: 
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c2U">) 


M-x 


2  e 


Cx(jo))    |(l  +— )2  +  sin2aj 


exp 


-l/sln_a  \ 
1  j.  x   j 


(131) 


Equation  (131)  is  then  transformed  Into  amplitude  ratio  and  phase 
shift  which  give  equation  (3k)    of  the  text. 
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The  general  extent  of  axial  dispersion  and  its  effect  on  the 
performance  of  isothermal  tubular  flow  reactors  are  investigated 
by  use  of  the  diffusion -type  model  and  four  sets  of  following 
boundary  conditions: 


B.C. (I) 


V- 


uCi  =  uCz_^0+  - 


S0z=L  =  ° 


3(32.) 


B.C. (II) 


uOx  =  ^Cz_0+ 


.(ML,-0 


B.C. (Ill) 


uCX  =  uCz^0+ 


lim  C(z,t)  =  0 
z-*-a> 


B.C 


uC]_  =  uCz_^o+ 


3(*2\ 


.(iv)  i 


lim  C(z,t)  =  0 

Z-+-00 


The  differences  and  the  relationship  among  the  steady-state 
concentration  profiles  and  conversion  and  unsteady-state  tran- 
sient response  curves  obtained  by  use  of  the  different  sets  of 
boundary  conditions  are  critically  re-examined  and  numerically 
evaluated.   Only  the  values  obtained  by  use  of  B.C. (I)  are  found 
to  be  consistent  with  the  limits  of  maximum  and  minimum  conver- 
sions, i.e.  conversion  in  a  plug-flow  reactor  and  in  a  complete- 


ly  stirred  tank  reactor  calculated  from  a  one -dimensional  trans- 
port equation0   The  solutions  obtained  by  use  of  the  other  three 
sets  of  boundary  conditions  are  found  to  be  inconsistent  with 
such  limits o   Hence,  unless  their  limits  are  recognizeds  the  use 
of  Bo  Co  (II)  j,  (III)  and  (IV)  should  be  avoided  0   Due  to  the  simple 
analytical  solutions  of  equations  in  exponential  form,  the  use 
of  BcCo(III)  and  (IV) ,  however,  offers  a  convenient  means  for 
computing  final  conversions  if  their  limits  are  recognized0 

An  illustration  is  presented  of  the  method  of  determining 
mean  dispersion  coefficients  from  dynamic  experiments 0   It  is 
recommended  that  the  calculated  values  should  be  compared  with 
the  experimental  data  to  determine  the  feasibility  of  predicted 
models  and  boundary  conditions „ 

The  response  curves  based  on  different  types  of  transient 
inputs  are  presented  for  all  four  different  sets  of  boundary 
conditions  both  with  and  without  chemical  reaction.   These  curves 
should  be  useful  for  practicing  process  control  engineers  and  for 
those  who  wish  to  extract  the  transport  parameter  from  dynamic 
response  experiments » 
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